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ABSTRACT 


Methods for transforming partial differential equations into forms more 
suitable for analysis and solution are investigated. The idea of a Gen- 
eralized Similarity Analysis is introduced and results applied to the equa- 
tions of boundary-layer flow. A thorough presentation of the application 
of continuous transformation groups to the problem of similarity analyses 
(reduction of the number of independent variables of an equation) is given 
with new and very general methods evolved for determining the nature of trans- 
formations . 



CHAPTER 1 


GENERAL CONCEPTS AND SCOPE OF THE INVESTIGATION 


1,0 INTRODUCTION 

The investigations to be presented in this report are an outgrowth of a 
continuing study of ways in which partial differential equations associated 
with problems of physical interest may be simplified through transformation 
of independent and dependent variables. The original motivation for this study 
gpew out of an interest on the part of investigators in the National Aero- 
nautics and Space Agency as to how rather general methods might be formulated 
for reducing the number of independent variables of a partial differential 
equation (similarity analyses ) . 

As the investigation proceeded, certain facts became clear. ' First it was 
noted that the usual types of transformation for simplifying partial . differen- 
tial equations of physical problems were usually of a rather special class. 

This led to the question of possibly generalizing the class of transformations 
noted. As will be seen significant success was achieved in this inquiry. 

t 

A second main area of investigation centered on finding methods of ana- 
lyzing very general types of partial differential equations other than spe- 
cific types such as viscous flow equations, etc. It was soon discovered that 
little of general nature could be uncovered in connection with the more com- 
monly used methods found in current literature. Most of these methods treated 
problems for which not only the differential equations but the boundary values 
were also specified. Transformations evolved were based on both the equation 
and the nature of the boundary values. There was one exception to this, how- 
ever, that gave promise of providing a very generalized method of analysis. 

This was the theory of continuous transformation groups. The method was by 
no means new. In fact, the basic ideas date back to the last century and are 
found in the work of the mathematician Sophus Lie. Moreover, the theory of 
groups has been employed quite extensively in recent times by investigators 
in the field of similarity analysis. Nevertheless, there was not in the lit- 
erature a truly in depth study of the applications of Lie’s theory to the 
similarity analysis of partial differential equations. Previous applications 
were quite limited in scope and it therefore seemed reasonable to find how 
far Lie’s ideas might be pursued in formulating a very general approach to 
similarity analyses. The results that have been obtained and that are pre- 
sented here are most encouraging. Very general methods of analysis have been 
formulated and the way made clear for further work. 
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The approach to he used in this report is the following. For the re- 
mainder of this chapter, a presentation of general concepts will he made, fol- 
lowed in Chapter 2 hy a discussion of the form of general transformations for 
partial differential equations of problems having a physical basis. The use 
of transformation groups for similarity analyses will be introduced in Chapter 
3, and basic ideas and applications presented. In Chapter 4, the methods em- 
ployed In Chapter 3 will be generalized and a set of rather powerful tech- 
niques for analyzing partial differential equations will be developed. 

1 . 1 GENERAL CONCEPTS 

The concept of similarity solutions is very old and dates back to such 
well-known workers as Buckingham^ and Bridgeman.^ Originally, the term sim- 
ilarity analysis implied a procedure for finding some Information about the 
solution of particular physical problem usually short of a complete mathemat-. 
ical answer. More recently, in the work of such analyists as Birkhoff^ and 
Sedov, ^ similarity solutions have incorporated rigorous mathematical techniques 
which result in solutions in the engineering sense, that is, numerical answers. 
Of particular importance, the more recent work has yielded solutions to non- 
linear partial differential equations which have been intractable to more 
standard solution techniques . In fact, exact solutions of the boundary-layer, 
equations of fluid mechanics are almost universally based on similarity methods. 

The various similarity solution methods can be broadly characterized as 
techniques which employ transformations of variables or parameters, or both. 
Typically transformations may linearize a problem (for example the Kirchhoff ' 
and ho do graph transformations), reduce partial differential equations to 
ordinary differential equations (for example the Blasius similarity trans- 
formations), transform the system to one already solved (such as the Mangier 
transformation), or perform some other reduction of mathematical complexity. 
Sedov 4 formulates a mathematical theory of similarity on the basis that "two 
phenomena are similar, if the assigned characteristics of the other by a simple 
conversion, which is analogous to the transformation from one system of units 
of measurement to another/" Thus, a certain analogy exists between the theories 
of dimensional analysis as formulated by Buckingham and Bridgeman and the geo- 
metric theory of invariants relative to a transformation of variables, a fun- 
damental : theory in modern mathematics and physics. Some authors have even 
found it convenient to define classes or types of similarity. Kline5 defines 
"external similitude" as the similarity between problems of given class (a 
transformation In terms of parameters alone) and "internal similitude" as a 
similarity characterized by a mathematical relationship between points inside 
a single system of a given class (for example, a transformation which reduces 
the number of independent variables). 

Inherent in any general theory of similarity, however, should be the rec- 
ognition of specifying, in some sense or other, the complete physical problem 
to be analyzed. In fact, the term "similarity" Implies a comparison between 
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two or more complete and recognizable systems. 

Thus, it is of prime importance to formulate a "complete" mathematical 
description of the physical problem under consideration before proceeding 
with a similarity solution. Gukhman^ characterizes a complete problem as one 
which is expressed in terms of the governing equations together with a body 
of information termed the "conditions for uniqueness of solution" which may 
contain not only boundary and initial conditions in the usual sense, but also 
possible auxiliary physical considerations such as conservation requirements. 

To be sure, it is not always possible' to prescribe the uniqueness conditions 
a priori fbr all problems that engineers are called on to consider. In fact, 
it appears that much of the work that has been done in the past to develop 
similarity analyses has been motivated by an initial statement of complete- 
ness, or lack of it. For example, the early work of Buckingham and Bridgeman 
on dimensional analysis require only a recognition of the pertinent variables 
which apply, without any statement being made concerning the governing equa- 
tions. Later workers such as Morgan, 7 Hansen,® Krzywoblocki,9 and Wecker and 
Hayes.10 investigated similarity methods by considering the governing equations 
first, and only examining the boundary and initial conditions as a later step, 
if at all. Another group of workers developed similarity methods by starting 
with a complete mathematical formulation and, thus motivated, to examine less 
complete (and more general) problems; see for example Coles, Abbott and 
Kline, and Gukhman. 

An examination of these earlier works show that the initial problem state- 
ment, as far as assumed completeness, determined to a large extent the kind of 
mathematical approach employed. The more information that was known, the more 
direct was the method developed for finding a similarity solution and, at the 
same time, the less general were both the methods and the conclusions (as re- 
gards "general solutions"). It is not suggested that this dichotomy is neces- 
sarily bad. The more general techniques, such as group theory methods, have 
produced powerful theorems and yield results with a minimum Df mathematical 
busy-work. On the other hand, the group theory methods are difficult for the 
average engineer to follow because their motivation is mathematical, not phys- 
ical, and this has inhibited their wide use. Also, somewhat amazingly, the 
more powerful mathematical techniques have been to a degree more restrictive 
in some of their aspects (such as the "class of assumed transformations") than 
the less elegant methods. 

By recognizing the differences in past motivation, and the resulting ad- 
vantages and weaknesses, the suggestion naturally arises that perhaps a 
marriage of the two approaches might be fruitful; this on one goal of the 
present investigation. In the following chapters, the two viewpoints will be 
examined in more detail than has been done in the past. In the following 
chapter a technique based strongly on physical reasoning will be used to 
evolve a very broad definition of similarity. It will be shown that by intro- 
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ducing the concept of generalized similarity,* a unified method can he de- 
veloped for deriving the majority of the -well-known transformations employed 
in fluid dynamics. 


*The term "generalized similarity" was found desirable because of the popular 
use of the term similarity to imply simply a reduction in the number of vari- 
ables of a given problem. 
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CHAPTER 2 


GENERAL TRANSFORMATION METHODS IN FLOW THEORY 


2.0 GENERAL CONCEPTS 

An appropriate change of variables is probably one of the most useful 
methods available for solving the partial differential equations of mathematical 
physics. The most general criterion for a transformation is simply to change 
a given problem into a simpler problem is some sense or other; that is, either 
to a form which will yield to more standard solution techniques, or possible 
to a form which has been previously solved in connection with a related or 
similar problem. Holograph transformations and conformal transformations are 
well-known methods for transforming given problems into forms which yield to 
classical techniques (such as separation of variables). Ames ^3 classifies the 
various transformation techniques into three . groups : transformations only of 

the dependent variables, transformations only of the independent variables, 
and mixed transformations of both independent and dependent variables. How- 
ever, all three groups have a common goal: to find a relation, or more spe- 

cifically, a basis of comparison, between different physical (or mathematical) 
problems. This broad concept of comparison is the definition of the term "gen- 
eralized similarity. " 

Generalized similarity might be applied in fluid mechanics to attempt to 
answer such questions as "Is there any similarity (i.e., basis of comparison) 
between compressible and incompressible flow problems, axisymmetric and planar 
flows, or in general, any more complicated and a less complicated flow?" By 
contrast, the usual term "similarity" as used by Birkhoff,^ Morgan,? Hansen,® 
Abbott and Kline., ^2 and others is defined in terms of independent variables 
of a problem.. Thus it may ultimately refer to a physical similarity within 
a given problem, such as similarity of velocity or temperature profiles. 

Much of the previous work on similarity was motivated by the desire to 
develop simple methods for reducing the number of independent variables. 

Out of this previous work came the realization that each of the proposed 
methods was based on an assumed class of transformations and the recognition 
that more general classes of transformations might lead to the solution of 
a wider class of problems. The first part of this chapter is concerned with 
seeking the most general class of transformations for particular types of 
problems. The two viewpoints ± discussed in Chapter 1 will be examined. 

First, the mathematical theory of transformations will be reviewed and it 
will be shown possible to postulate the so-called primitive transformation 
as the most general form of a class of transformations (under a certain as- 
sumption). Second, a separate approach, based on postulating the "complete 
physical problem, " is examined for the special case of laminar boundary-layer 
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flows and it is shown that this approach also yields the primitive transforma 
tion. 


The second part of the chapter deals with the development of a technique 
for employing the primitive transformation to find the form of the variables 
for a wide variety of generalized similarity problems. 


2.1 TYPICAL TRANSFORMATIONS USED IN FLUID MECHANICS 

As a means of developing some insight into the question of "the most gen 
eral class of similarity transformations" and the concept of "generalized 
similarity, " Table 1 was compiles of examples of as many different types of 
known transformations as could be found to represent the field of fluid me- 
chanics. 


TABLE 1 


Trans format ion^ 

Transformed Independent 

Variables 

Similarity 

i(x) 

= x, T)(x,y) = y y(x) 


Meksyn-Gortler 

l(x) 

X 

= / u x (x)dx, T](x,y) 

o 

= u 1 (x)y/ n/x 

von Mises 

l(x) 

= x, Ti(x,y) = t(x,y) 


Crocco 

t(x) 

= x, T](x,y) = u(x,y) 


Mangier 

S (x) 

= / X r§(x)dx, Tj(x,y) 
o 

= r 0 (x-)y 

Stewartson 

!(x) 

X 

= / a 1 (x)p 1 (x)dx, 



o 


X 

Tl(x,y) = / a 3 .(x)p(x,y)dy 
o 

Dorodnitsyn |(x) = / Pi(x)dx, n(x,y) ■= / p(x,y)dy 

o o 

Similarity Rules of High Speed Flow (see, for example Ref. 16): 

(i) Prandtl-Glauett 

(ii) Gothert |(x) = x, q(y) = By 

(iii) von K£rrn£n Transonic 
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By comparing the various transformations in the table, an interesting 
conclusion can be made: all of the transformed independent variables are of 

the general form | = f(x) and T] = T](x,y) (that is, one of the new variables 
is a function of only one of the original variables). This realization 
raises the question "Is the general transformation |(x), Tj(x,y) the most gen- 
eral class of transformations that need be considered for physical problems?" 
In an attempt to answer this question, it will first prove useful to review 
the mathematical theory of transformations . 


2,2 THE PRIMITIVE TRANSFORMATION 

The following theorem can be found in the mathematical literature of gen- 
eral transformation theory (see, for example, Courant-^ ) : 

THEOREM: An arbitrary one-to-one continuously diffe'rentiable transformation 


i = l(x,y), t) = T)(x,y) 


of a region R in the x,y-plane onto a region R 1 in the £,t} - plane can be re- 
solved in the neighborhood of any point interior to R into one or more con- 
tinuously differentiable "primitive" transformations, provided that through- 
out the region R, the Jacobian 


J(£,T]) 


d(g,T)) 

d(x,y) 


Ix’ly ~ iy'Hx 


differs from zero. A "primitive" transformation is of the form 


| = g(x), tj = ri(x,y) • 


Now one-to-one transformations have an important interpretation and ap- 
plication in the representation of deformation or motions of continuously dis- 
tributed systems, such as fluids. For example, if a fluid is spread out at a 
given time over a region Rand then is deformed by motion, the motion of the 
fluid is described by the coordinates in the physical R-plane. If the fluid 
motion in R is characterized by the coordinates x,y, then the corresponding 
motion in the transformed region R 1 is characterized by coordinates §,t}. The 
one-to-one character of the transformation obtained by bringing every point 
x,y into correspondence with a single point g,T| is simply the mathematical ex- 
pression of the physically obvious fact that the fluid motion in the physical 
R-plane must remain recognizable after transformation to the transformed R'- 
plane, i.e., that the corresponding motions remain distinguishable. 
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Physically, since most "practical" transformations of interest for solving 
physical problems should he one-to-one,* that is, have a unique inverse (ex- 
cept possibly at a finite number of singular points), it appears the primitive 
transformation or the resolution into primitive transformations, should be con- 
sidered in the search for "the most general class of transformations." 

For the sake of completeness, the following two properties of the primi- 
tive transformation are noted: 

(i) If the primitive transformation 


| = |(x), T] = T](x,y) 


is continuously differentiable, and its Jacobian 




d(g.»Tl) 

d(x,y) 




’^X 


differs from zero at a point P(x 0 ,y 0 ), then in the neighborhpod of 
P the transformation has a unique inverse and this inverse is also 
a primitive transformation of the same type. ! : 

(ii) For primitive transformations ^ the sense of rotation in the x^y- 
plane is preserved or reversed in the -plane according as the 
sign of the Jacobian is positive or negative 3 respectively . 

In summary^ the primitive transformation appears to he the most general 
class of transformations that need he considered^ purely* from a mathematical 
viewpoint^ as long as it is required that a unique inverse of the transforma- 
tion must exist. In the next section,, the question of the most general trans- 
formation will he approached from a fundamentally physical viewpoint and it 
will he shown that again the primitive transformation appears to he a require- 
ment from conclusions based on uniqueness arguments. 


2.3 REQUIREMENTS IMPOSED BY THE PHYSICAL DESCRIPTION 

The question of the most general class of transformations will now he 


*An exception to this rule; for example^ is the von Mises transformation of 
Table 1, which is singular along the x-axis. However ^ this transformation 
is used for computational (not physical) reasons^ and hence its use is 
motivated by a completely different line of reasoning than that under con- 
sideration here. 
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examined from a more physical veiwpoint. Some mathematics -will he involved, 
of course, hut the physical problem will he kept near at hand and frequent 
reference to it will he made throughout the analysis. 

Because in what follows, hy definition, involves a particular physical 
problem, it will he convenient to introduce such a problem. As an example, 
consider the following equations describing the two-dimensional flow of a 
laminar incompressible boundary layer: 


( 2 . 1 ) 

dp 

dx 


or an equivalent single equation in terms of the stream function; 

St _ chjr Sft _ Sjt _ dp (2.2) 

Sy SxSy Sx Sy 2 Sy 3 dx 

For convience in the development, it is assumed that all of the variables are 
nondimensional. The objective is now to consider very general transformation 
of variables and see what conditions must be met by this transformation so as 
not to violate any known physical properties of the problem. 

Let us begin by specifying a transformation of variables. For the inde- 
pendent variables, no restrictions at present will be placed on the assumed 
form; thus the form is specified simply as 


du . bv n 

TT + = U 

ox oy 


y 


Bu 


U Bx + v dy 


£u a 2 u 


c>y 2 


I = l(x,y), q = T)(x,y) . 


(2-3) 


The dependent variable ljr (the stream function) is also transformed to a new 
function, say ¥. Previous work on fluid flow tranformations has often assumed 
that the two stream functions, ¥ and ¥, should be the same at corresponding 
points and hence that streamlines in one plane are tranformed into streamlines 
in the other. However, this restriction is found to be unwarranted in many 
problems and will be avoided here. Instead, the relationship between i and ¥ 
will be specified in some weak sense by the form 


¥(|,ti) = g(x,y)t(x,y) 


(2.4) 
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It should, he recognized that this assumption is not trivial and other forms 
could be considered. Nevertheless* Eq.. (2.4) represents a more general case 
than usually employed and attention will he restricted to this form in the 
present report. 

The next step is to carry- out the transformation of Eqs. (2.3) and (2.4) 
on the left-hand side of Eq.. (2.2). Application of the transformations (2.3) 
and (2.4) to the boundary -layer equation is not new and was considered pre- 
viously by Colesl5 and others. However* the generalized similarity technique 
was not under consideration by the previous authors . 

The results of the transformation of Eq.. (2.2) is given in Table 2. For 
convenience in interpreting the physical terms after transformation* trans- 
formed velocities U*V have been defined by 


U 



d£ 


The transformed side of Eq.. (2.2) given in Table 2 is essentially un- 
manageable in its present form. While it is unclear at the present time what 
conditions are either necessary or sufficient to ensure any particular math- 
ematical behavior* it is interesting to consider an argument proposed by 
Coles. -*-5 Coles’ argument is based on the a priori requirement that the trans- 
formed flow outside of the shear layer (i.e.* for large values of y) should 
conform both physically and formally to the original flow. Thus* he argues, 
since the, physical flow is bounded for large y and is* in fact* at most a 
function of x* then the transformed side of Eq.. (2.2) must behave in a sim- 
ilar fashion. His argument is then that the substantial derivative terms 
and the U 2 terms become at most functions of £ for large y* whereas the re- 
maining terms behave either like y or y 2 . His conclusion is to require the 
r] and rj 2 terms to vanish identically* which can be accomplished by requiring 
that g = g(x)* | = |(x)* and r)y = y(x). At the present time* about all that 
can be said is that this assumed form is one possible form of a transforma- 
tion which will preserve a certain sense of physical correspondence between 
the given and transformed flow; this form of the transformation will be em- 
ployed throughout the remainder of the present chapter. 

In summary* the postulated general transformation* based on the physical 
mode* is : 


¥(£*t]) = g(x)!(x*y) 

£ = i(x)* n = y(y(x) ) 


(2.5) 
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TABLE 2 


Transformation: = g(x J y)t(x > y) ; g = g(x,y),t| = rj(x,y), U = ^ } V 

dr) 


dp = st s 2 t _ Sj/ sft _ sfi (r Su su\ / 

dx Sy SxSy Sx Sy 2 Sy 3 V SrJ/ y 

■ 

• *(M, • ! &X b - 3 ■»] 


y»^ 


+ u2 jUj/gj) . OT J(gy/e.n)w(y/s,i) + ^ 


J(6v/g,l) 


M « -2 

S| LU /y 


g yyy + 2 ^ . y r l z ) + ^lz) i J 


+ fiz') * .fzS+^lzV 

S| IA s Jy y g^ \ g ,/y_ 


Su A 


Z . _5L^J'2L 


stiLU y y n y " g 2 V + U )y\ 

ifu ril n + 2 ffia TI + li 

S| 2 l_g y g J g 

sfu if + sfv ^y - + 2 ^y^y 

Sri 2 s Sn 2 Lg y g J 


i T |7^ U n ^ V J (g> T l) +( n - 5v ? ) J (g>l). 4-f !“Z 

LUV ^ VSg v s thj g3 + ^s. 




11 



It is of interest to note that the resulting form of the independent variables 
is a primitive transformation as discussed in the last section, hut it has 
"been shown by considering a particular problem and relying on physical argu- 
ments that it may be possible to restrict the form of T] to be a linear func- 
tion of y for this particular class of problems (boundary-layer flows). Thus, 
the two approaches, mathematical and physical, have lead to the same con- 
clusion regarding a postulated "most general class of transformation." Of 
course it is necessary, as previously pointed out, to emphasize that there 
may be some cases of practical interest which will lie outside the realm of 
the mathematical conclusions presented here. However, it has been found that 
all of the cases listed in Table 1 are satisfied by the transformation of Eq. 
(2-5) (with the exception of the von Mises transformation previously discussed). 

It is now worthwhile to return to Table 1 for a moment. Recall that from 
Table 1, it was noted that a wide class of different transformations were all 
primitive transformations. Further, note that one of these transformations, 
the similarity transformation, has been shown to yield to simple general anal- 
ysis for its derivation for particular problems. In fact, at the present time 
there are two types of similarity analyses that are founded primarily on sim- 
ple transformation theory: the free parameter method of Hansen® and the 

separation of variables method of Abbott and Kline. ^ It is thus of interest 
to speculate on the following question: Would it be possible to derive all 

of these "generalized similarity" transformations by the same technique that 
is used to derive the similarity transformation? The Generalized Similarity 
Analysis was developed as an attempt to answer this question. 

This new method is based on a single assumption concerning the admissible 
class of transformations of the independent variables, namely that the assumed 
transformation should be one-to-one (that is, have a unique inverse). On the 
previous pages it was shown that this requirement would lead to the primitive 
transformation, and further that a particular primitive transformation may be 
obtained from physical arguments for a particular class of problems, namely 
i = |(x), T) = yy(x) for the boundary -layer flows. Of course the generality 
of the present ideas goes beyond a particular case, such as boundary- layer 
flow analysis and should be appliable to a wide class of problems. 


2.4 THE GENERALIZED SIMILARITY ANALYSIS 

The development of the generalized similarity analysis was motivated as 
an extension of the method of finding similarity variables (i.e., the reduc- 
tion of the number of independent variables) to the problem of finding a trans- 
formation of variables which will convert a given physical problem into an 
alternate problem under certain prescribed conditions. For example, the pre- 
scribed condition for a similarity solution is that the number of independent 
variables must be reduced. By contrast, the prescribed condition for the 
Mangier transformation is that the axisymmetric boundary-layer equations are 
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to "be transformed into the planar 'form of the equations, and so forth, with 
similar statements for the rest of the examples in Table 1. 

There are three distinct. ..steps to the generalized similarity analysis. 
These steps are: 

(a) The general mathemathical theory of transformations states that any 
continuous one-to-one transformation can he resolved into one or more prim- 
itive transformations of the form 


| = !(x), T) = T)(x,y). 


Thus, a primitive transformation form is assumed a priori, where it is recog- 
nized that the general-' analysis has the possibility of being repeated more than 
once, depending on the particular problem at hand. 

(b) The given equation, transformed under a primitive transformation to 
the new independent variables (j-,T]) is required to satisfy the state require- 
ment; for example, that for a similarity analysis, the transformed equation 
should be a function only one of the new variables. 

(c) Simultaneously the boundary conditions for the given equation are 
required to be satisfied when expressed in terms of the transformed variables. 

' These three steps will be shown to completely and uniquely determine the 
explicit form of the 'new variables £ (x) and t] (x,y ) if, in fact, the original 
problem and associated boundary conditions do admit a generalized similarity 
solution. 

1 As an example of the method, two problems will be examined in detail. 

First the classical Blasius flat plate boundary-layer problem will be solved 
to give a relatively simple motivation of the basic ideas. Second and more 
difficult Mangier transformation will be derived as an example of the broad 
applicability of the method. 


Example 1 


The boundary-layer equations for a steady, two-dimensional laminar flow 
with a zero pressure gradient (Blasius flat plate flow) can be written in the 
following form in terms' of the stream function: 


chjr cfy 6 2 ijr _ 5 3 t|t 

6y dxdy 6x £y 2 £y3 


(2.6) 
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(all variables are dimensional, where v is the kinematic viscosity) with the 
boundary conditions 


y - 0: |t . 0, ft = 0 

dy Sx (2.7) 

y oo : ^ 4 . -> constant = u 0 . 

oy 

To make the problem statement complete, it is necessary to state the given re- 
quirement for the transformation x,y |,T): 

Requirement: under the transformation, the given equation and boundary condi- 

tions, (2.6) and (2.7)> are to reduce to a function of a single 
independent variable (i.e , the similarity variable). 

The three steps of the generalized similarity analysis are carried out in order 
as follows: 

(a) Assume a primitive transformation 


| = g(x) and T) = T}(x,y). 


Since a boundary-layer problem is under consideration, it was shown in the 
Section 2.3 that it is sufficient to choose the more specific transformation 
for T) as T) = yy(x), and further, for the present case of a similarity solu- 
tion, | can be assumed in the simple form £ = x without loss of generality 
because the final result, by definition, is to be a function of the single 
variable T](x,y), and not both | and T).* For the dependent variable f(x,y), 
the transformed variable ¥(|,r|) is assumed in the form 


¥(£,r| ) = g(x)\|f(x,y) 

as discussed in Section 2.3. Thus, the assumed transformation for the depen- 
dent and independent variables becomes 


*The question naturally arises, why not assume the form % - y, T) = T|(x,y)? 
This alternative is easily eliminated for the present case by carrying out 
the corresponding analyses in an analogous manner as that presented here; 
it is then found that this assumed form will not satisfy the given boundary 
conditions . 



T, g(x)t(x,y) ■ 
i = X, TJ = yy(x) • 


( 2 . 8 ) 


Performing the transformation (2.8) on the given Eq. (2.6) yields the 
following results: 


8jjr _ 8t fy 
8y . ~ 8ti g 


8ft 

by 2 


8ft 

8ti 2 g 


8 3 t 

83^ 


8ft fi 

8t) 3 g 


8t 

8x 


^g^_ + 8tl + 8tfx 

g 2 8| g - 8r| g 


8 2 t = 8tAxy g» Y _8ft fy 8 2 t 'Hy'Hx 

8x8y 8ri \ g g 2 y g + ^2 g 

where primes are used to indicate total derivatives. Substituting these re- 
sults into Eq.. (2.6) yields: 


8 s ¥ ' /g 1 1 . _ + At 8t 8 2 t \ 1 _ f 8ft gh_l 

^ 3 \g T) Al Y 2 ^T\ 8|8r| / gTjy 8r| 2 g T]y 

■with boundary conditions : 

tj » T|(x,0) * O: ~ n y (x,0) = 0; U^ 0 ) = 0 

rj = Tj(x^y) 00 : ^ — y . fey * 00 )_ constant = u D • 

8r] g 


(■b) Imposing the stated condition of similarity by requiring that 

t = t( T| ) 
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that is, that the transformation must reduce the number of variables, and 
noting that r\y = y(x), the transformed equation becomes: 


Vf 1 ' 


+ (-&L. 

\^7 


- lL 


r 2 g 




g T 
2 
g 


= 0 


with boundary conditions 


tj = 0: T'7 =0, -¥ + ¥'7 = 0 

g 


rj ■+■ oo : ¥' — ■> u 0 • 


( 2 . 9 ) 


( 2 . 10 ) 


(c) So far the functions g(x) and y(x) are unspecified. However, only 
particular forms of these functions will satisfy "both the boundary conditions 
and the stated requirement that Eq_. (2.9) must be a function only of T) . Ex- 
amining the last boundary condition (2.9) yields 


f ' (00) iM constant = Uq . 
g(x) 

The only way that this can be true is if 


7 ( x ) = u 0 g(x) . 


Also, the only way that Eq.. (2,9) can be a function of the single variable rj 
is if the coefficients of ¥ and its derivatives are constant. Hence, letting 


2 -2 — = constant = c x 

g 2 7 u 0 g 3 

then integration of this ordinary differential equation yields 


g 


[- 2 c 1 u 0 (x+x 0 )]" 1//2 
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where Xq is an undeterminable constant of integration,* and letting the arbi- 
trary constant c x take the value c x = -v/2 yields 


g 


[ vu 0 (x + 


Xo)3 


- 1/2 


7 


Uq [ VUq ( X + Xq)] 


- 1/2 


Evaluating the coefficient of ¥' 2 : 



Ci and 


7' 
7 s g 


ttQS 1 = c 


Thus this coefficient is identically zero** and the final form of similarity 
solution is: 


Y"i + i. y?" = 0 

2 

¥(0) = ¥'(0) = 0, y«(oo) 1 . 


( 2 . 11 ) 


In summary, it has been shown that a unique choice for the transformation 
variables can be found for the similarity solution of the Blasius problem by 
carrying out the three steps of the generalized similarity analysis. The next 
example will show that the method can be applied to a much more- generalized 
problem, that of finding the conditions for which an axisymmetric boundary 
layer is similar to a planar boundary layer. 

*The meaning of Xq is clear from the physips of the problem] its magnitude 
is. determined by conditions at the leading Sdge of the plate and cannot be 
derived from boundary -layer theory alone. Thus the appearance of Xq is a 
sufficient proof that the Blasius solution is a correct asymptotic down- 
stream solution for laminar flow over a plate which holds independent of 
the initial conditions specifying the plate leading edge. 

**The coefficient of f^ 2 in Bq.- (2.9) is not independent of the relation 
g , /g 2 7 - c x . This can be seen as follows: for Simplicity, let 7 = l/6 and 
g = l/h. ThSn g’/g 2 7 = -h*8 = c x and letting (g’/g 2 ! - ! '/7 2 g) = (-h'5+S’h) 

= C2, then S'h is a constant, say c 3 . The two equations h'5 = -c x and 
5’h = c 3 maybe combined to yield 55" = (c x /c 3 )S 12 = 0 which has the solu- 
tion 5 = |ax + b | 1 /(1 _g i/ c 3) f 0r Cl ^ C3 an< 3_ 5 = ae^ x for c x = c 3 . This 
provides a general solution, however we cannot evalute c x or c 3 without 
another condition. The conclusion, as found above, remains the same; the 
required condition must come from the boundary conditions . 
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Example 2 


The well-known Mangier transformation D serves as a good example of the 
broad meaning of the idea of generalized similarity because the transformation 
answers the question "Under what conditions is an axisymmetric boundary-layer 
flow similar to a two-dimensional planar flow?" The answer lies in finding 
a transformation between the variables describing the two types of flow. The 
generalized similarity analysis is formulated in the same way as the simple 
similarity analysis of the preceding example; that is, by specifying the equa- 
tions (and boundary conditions) and the stated requirement to be satisfied 
by the transformation. 

The governing equations for a thin laminar axisymmetric boundary layer 

are: 


d(g 0 u-) d(r 0 v) 
dx + Sy 


(2.12) 


u ^ + v ^ 
ox dy 


U 1 


dui + v d u 
dx by 2 


(2.13) 


where r Q = r D (x) is a given quantity (r D specifies the body shape relative to 
the axis of symmetry) and v, the kinematic viscosity, is a constant. The gen- 
eral boundary conditions are: 


y = 0 : u, v = 0 

y -> oo : u + u x (x) • 


(2.14) 


The problem statement is completed by writing down the required mathematical 
form of the desired equations, namely: 

Requirement: under the transformation u,v ->• U,V and x,y -* |,T) the trans- 

formed equations are to be in the form of a planar boundary- 
layer flow; that is, of the form 


8u &V 

+ 


0 


(2.15) 


u ! 


+ v 


8u 

St) 


= Un 


dU-i 

d| 


+ V 


aPa 

dr ] 2 


(2.16) 
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The three parts of the analysis are as follows: 

(a) Assume a primitive transformation. Again, since the problem involves 
the boundary- layer equations, it is sufficient to assume the same general 
form of the transformation as given fn Section. 2.3, that is 


*(t,Tl) = g(x)\|r(x,y) 

I = I (x), n = yr( x ) 


(2.17) 


where, for the present case, the stream function f (x,y) is defined by the equa- 
tions 


1 dvlr 

u = — , 

r Q dy 


v = 


1_ d£ 

r 0 Sx 


so as to satisfy Eq.. (2.12). For the transformed flow to be in the form of 
planar equations, the transformed stream function ¥(t, T l) is defined by 


u = 5*1 , v = - ^ 


dr] 


*6 


Using these definitions and the transformation (2.17), the velocities and 
their derivatives become 


v(x,y) = 


u(x,y) = 

7 

u(i,n) 


r Q g 


*(l/l) 

+ l25_ 

v(i,n) - 

g 

r oS 

r o§ 

du _ 

_2_ - 


Sy 

r 0 S 

dri 

dfu 

JL- y 

2 d 2 U 


r o§ 

drj 2 


dU _ ( 7 V TT , _y_/dU dU \ 

dx \r Q g/ r Q g \ d| dr) 
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Substituting these expressions into Eqs. (2.12) and (2.13) yields, after re- 
arrangement : 


J(t,n) 


av\ 

V 


0 


lx 



( 7 r o g ) v 


+ 



(uf-u 2 ) _ £l f ^ . 

g di) 


(2.18) 


(2.19) 


It is found that the continuity Eq. (2.l8) is invariant under the transforma- 
tion (2.17) as long as the Jacobian J(g,q) = y| x differs from zero. Thus Eq. 
(2.l8) does not force any requirements on the transformation. However, in 
comparing Eqs. (2.16) and (2.I9), it is seen that for the latter equation to 
fulfill the given requirement of being identical to the planar form given by 
Eq. (2.16), it is necessary that 


lx = 7 r oS 


(2.20) 


and 



A sufficient condition for Eq. (2.21) to be satisfied is for 


( 2 . 21 ) 


and 



constant = c x 


g = constant = C2 • 


( 2 . 22 ) 


(2.23) 


Although Eqs. (2.22) and (2.23) may not be necessary conditions for the satis- 
faction of (2.21), they at least provide one solution for the given require- 


20 





merit and this is usually satisfactory from an engineering veiwpoint. Since 
there are no further requirements to restrict a choice for the constants c x 
and cs, they are normally chosen for dimensional reasons to assume the values 
Ci = 1 and g = l/D where D is an arbitrary reference length. Hence, the final 
form for the transformation becomes: 


*(S,n) = ” t(x,y) 


(2.24) 


I 


r x r|(x) 

l ~W~ 


dx 


(2.25) 
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(2.26) 


(b) It has already been shown that for the present case, the unique form 
of the transformation, Eqs, (2.24), (2.25), and (2.26), is determined by re- 
quirements on the differential equation alone. Thus, the boundary conditions 
do not provide any additional information and, in fact, they are found to carry 
over directly as follows: 


T) = 0: U,V = 0 

Ti oo: U -*■ ti(i) 


(2.27) 


The analysis is thus complete, Eqs. (2.24), (2.25), and (2.26) being known as 
the Mangier transformation. 


2.5 CONCLUDING REMARKS 

It has been shown that the primitive transformation appears ot be the 
most general class of transformations necessary to provide a transformation 
with a unique inverse. This result was obtained from the general mathematical 
theory of transformations, but was also supported by a physical argument for 
such cases as boundary-layer problems which showed that the transformation 


I ( i , tj ) = g(x)\|r(x,y) 

% = l(x), r| = yy(x) 
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is possibly sufficient to ensure proper "behavior of the equations. 

The generalized similarity analysis was then introduced to solve a wide 
class of problems which could he formulated as a comparison between two given 
subproblems. Two examples were given, the Blasius similarity problem and a 
derivation of the Mangier transformation, however all of the cases given in 
Table 1 can be derived in s aimilar fashion. 

A few comments can be made concerning the role of the function g(x) ap- 
pearing in the transformation of the dependent variable for the boundary- layer 
equations. In certain cases, the value of g(x) will be unique ly det ermined 
by the problem. For example, g(x) is proportional to nTx and */{• (x) for the 
similarity and Meksyn-Gortler transformations, respectively, and g is found 
to be constant for the Mangier transformation. However, in some cases, the 
choice for g ia arbitrary; for example, g may take on any value for the von 
Mises transformation. Further, for the case of comparing compressible and in- 
compressible forms of the boundary-layer equations, different choices for g 
lead to different, but nevertheless useful, results: Stewartson chooses g = 

constant and Dorodnitsyn chooses g = V|y = y(x). Coles discusses the signif- 
icance of g for the difficult problem of compressible -incompressible trans- 
formations of the boundary-layer equations for turbulent flow in Ref. 11. 

In summary, the ideas presented in this chapter are based on, or more ap- 
propriately, motivated by a physical description of a problem which is in some 
sense complete. Depending on the particular case under consideration, com- 
pleteness may imply a knowledge of all necessary boundary conditions, or pos- 
sibly only a statement of a particular requirement of the transformation. In 
any case, a fairly specific problem formulation is implied. 

In the next chapter, a different technique will be examined which focuses 
attention on a more narrow application; the simple similarity problem (in the 
sense of reducing the number of independent variables). This technique, the 
group theory method, being less encumbered by statements of broad generality, 
will prove to yield very elegant and powerful mathematical results for find- 
ing similarity solutions for a wide range of applications. 
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CHAPTER 3 


SIMPLE GROUP THEORY APPLICATIONS 


3.0 INTRODUCTION 

In this chapter we will consider one of the more widely used methods of 
reducing the number of independent variables of a partial differential equa- 
tion, the group -theoretical method. The foundation of the method is contained 
in the general theories of continuous transformation group that were intro- 
duced and treated extensively by Lie in the latter part of the last century. 

G. Birkhoff^ was one of the first to apply this concept in searching for sim- 
ilarity solutions of partial differential equations . Subsequently, Morgan^ 
investigated quite thoroughly the mathematical theories involved and essen- 
tially completed the development of the method now generally referred to as 
the group -theoretic method. 5 

In the classical group-theoretic similarity analysis, a one -parameter 
linear or spiral group of transformation is employed.* The variance of a 
given partial differential equation under the transformation group and the in- 
troduction of new invariant variables, enables the number of independent var- 
iables to be reduced by one. Once the theories involved are understood, the 
process of obtaining similarity solutions is considered to be the simplest as 
compared with other techniques. Extension of this method to multi -parameter: 
linear and spiral groups of transformation has been developed in recent years, 
and will be discussed in this chapter along with the one-parameter approach. 


3-1 ONE-PARAMETER GROUP -THEORETIC ANALYSIS 

In this section the topic of group-theoretic analysis will be introduced 
by restricting attention to one -parameter methods. The material will have the 
twofold purpose of showing how the one-parameter method is used in similarity 
analyses and will provide a background for material in the remainder of this 
chapter and Chapter 4. By no means does the one-parameter similarity analysis 
plumb the depths of Lie's theory of continuous transformation groups. As 
will be seen, only two simple types of transformation groups will be con- 
sidered — the linear and spiral groups. In the next chapter, a far more gen- 
eral method is presented that permits the construction of a greater variety 
Of transformation groups. This section primarily reviews the work of Birk- 
hoff^ and Morgan. 5 


*See Ref. 3 for a general introduction to group theory and this method. 
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3.1.1 Absolute Invariants 

The first step in introducing the concept of one -parameter transformation 
groups is by establishing the group property of a transformation of coordinates . 
In general, if a transformation is defined by 


X 1 = f 1 ( x ' 




a', . . 


,n'i 


(3-1) 


then there exists (m-l) functionally independent absolute invariants 

(x' , . . . ,x m ) vhere j=l, . . . , (m-l ) . By definition, an absolute invariant has 
the property 


5J(x*,...,xF) = t ( x x m ) 


(3.2) 


Foilwing Morgan, 5 we restrict ourselves to the one-parameter group* 


T a : X" 1 = f^ - (x 1 , . . . ,x m ja ) 

yj = hJ(y', . . .,y n ;a) 


(3-3) 


where i=l, ... ,m (m>2 ) and j=l, . ..,n (n>l). 

For the subgroup X-*- = f^(x', . . .,x m ;a), there exists (m-l) functionally 
independent absolute invariants 


T)^ ( x 1 , • • . , x®' ) , k— 1, « . • , (m— 1 ) . 


(3-M 


For the group T a as a whole, there are (m+n-l) functionally independent abso- 
lute invariants. We therefore add the following to the list: 


g / g(x’, . . .^jy', . . .,y n ) , 
Jt=l, « u . > :n'v ’ v 


(3.5) 


*The variables x^ and yi will shortly be identified as independent and dependent 
variables. For now, they are just general variables. 


26 



3-1.2 Morgan's Theorem 


Assuming that yj are the dependent variables and x 1 the independent var- 
iables in the one-parameter group of transformations 


T a : X 1 =f i (x',...,x ffl ja) 

= h«3(y', . - -,y n ;a) , 


(3-6) 


the system of partial differential equations of order k 

Ac,..., = 0 

' &(*')* S(x m )V 


(3-7) 


is invariant under this group of .transformation, T a , if each of the ^ is "con- 
formally invariant"' under the transformation T k . This means that 

cA 


A ( vi vHl.yi y n . 8 7' 

r-il A ) • • ' j a > 1 > • ■ • > 1 > 7" ) • • • s 


, k v. a k Y n 


x d(X') k 

= f(x', . . .jX^y', - - .y n ; . . 


d(X m ) 


-m\k 


•vk n 

•n- . 2 I 


S(x m ) 


k J 


; a 


/ d k v n 


(3.8) 


In particular, if F = F(a) = 1, <f ^ is said to be "absolute invariant" under 
this group of transformations.* 

The above result leads to the very important theorem of Morgan: 5 

Theorem : Suppose that the forms ^ are conformally invariant under the group 

T a « Then the invariant solutions of =0 can be expressed in terns of solu- 
tions of a new system of partial differential equations 


/ ^f ^ k p 

Plij • • • > r lm-i;Fi, - • • ,F n j . . . ; — ^ ) = 0 

cTli 


(3-9) 


The r)j_ are the absolute invariants of the subgroup of transformations on the 
x 1 alone and the variables Fj are such that 


*See Ref. 3- 
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t> * 




(3.10) 


(Hi, • « • jTfrn-i) = function of (x'. 


An example is given in the next section to illustrate the method. 


3.1.3 Application of One-Parameter Group Theory to Boundary-Layer Equations 
Consider the steady, two-dimensional laminar houndary-layer equations: 


8u chr . 
dx by 


0 


(3.11) 


u Su + v 5u = udU+v^u 
bx by dx 3y2 

subject to the boundary conditions 


y = 0: u = v = 0 

y =. oo; u = U(x) 


( 3 - 12 ) 


Group theoretic! methods will now be employed to reduce 'Eqs* (3.1l) and 
(3.12) to ordinary differential equations* To this end, two groups of trans- 
formation (namely, the linear and the spiral groups) are considered. The 
method by Birkhoff and Morgan serves as the basis for the analysis that fol- 
lows . 


Case 1: The linear group. 

The linear group of transformation is defined as 


a«i- . Qh- .«3- 

x = A x, y-Ay, u = Au 


v = A tt4 v, U = A^U. 


(3.13) 


where are constants and A is the parameter of transformation. The 

function U will be considered as an independent variable with the defining re- 
lation U = U(x) . 
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Under this group of transformation, Eqs* (3*11) and (3*12) become 


and 


chi + a 4 -a 2 8v 
3x 8y 


(3.14) 


p2a s -a, x - ^ + A a 3 +a 4 -ai 2 ^ ^ 

Sx 8y 


A 2a 5 -a x 


U §1 + a «3- 2 «2 

dx 



(3-15) 


It is seen that the differential equations will he invariant before and 
after the transformation if the powers of A in each term are the same. Thus, 
we get 


a 3 - a x = a 4 - a 2 


(3.16) 


and 


2a 3 - a x = a 3 + a 4 - a 2 = 2 c% - a x = a 3 - 2a 2 . (3. IT) 

Solution of Eqs* (3«l6) and (J.l?) then g;ives 

a 3 = a 5 = a a - 2a 2 , a 4 = -a 2 , (3*18.) 

Knowing the relation among the a's, the absolute invariants can be ob- 
tained by eliminating the parameter of transformation A. Thus, noticing that 


II 

^ , a - (Xz/oti 

X 

U 

J. u 

x~zoc 

X 

- 1 - 20 : 

X 

Y 

V 

“ 

—-a 

X 

U 

U 

X 1 ~20£ 

- -i-aa ' 
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these combination of variables are readily shown to be invariant under the 
linear group of transformations and so are absolute invariants. According to 
Morgan's theorem, Eqs. (3. 11) and (3. 12) can now be expressed in terms of 
these invariants; or in other words, these invariants are possible similarity 
variables. Therefore, we put 

71 = x i-m/2 ’ = x® ’ = x m-x/2 (m=l-2a) (3»19) 

and 

hCil) = Ci = \ (3.20) 

X 


where f and g are functions of tj and h(Ti) = c x is a constant since U, the 
mainstream velocity, is a function of x only and thus cannot be a nonconstant 
function of T) which would introduce dependence on y. 

From Eq. (3*20), we get the form of mainstream velocity for the establish- 
ment of similarity solutions, i.e., 


U(x) = c^ . 


( 3 * 21 ) 


We may now make a check, the transformation of the boundary conditions;* 
It is seen that the boundary conditions are transformed to: 


r|=0: f s g = 0 


q = 00; 



Cl . 


Thus, for constant values of tj constant values of f and g result. It is now 
a simple matter to transform the differential equations. Thus, Eqs* (3*11 ) 
and (3.12) are transformed to 


mf - — qf • + g» = 0 

2 


( 3 * 22 ) 


and 


mf 2 _ IzB. rjff * + gf;« = me 2 + vf" 


(3*23) 
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The boundary condition are those previously stated. 


Case 2: The spiral group. 

A spiral group of transformation is defined hy 


x = x + p^b, 


y = 


- p 2 l3 

y e , 


- Ps* 
u = u e 


v = v e^ 4b , U = U e Psb . 


(3-24) 


Following the same steps as given in Case 1, we find that the parameters 
are related by 


&3 - &5 ~ ~&2 > P 4 - ~2p2 


(3-25) 


and, after elimination of b, the absolute invariants are: 


JL 


- 1 ex 
e 2 p 


f = 


u 

,P* 


1 

e2 


px 


(p=p 3 /pi: 


(3.26) 


and 

h(ri) = c 2 = . (3*27) 

eP x 

Again, c 2 must be a constant. Thus, the mainstream velocity should be of the 
form 


U(x) = C 2 e^ X 


(3-28) 


for similarity solutions to exist. Checking the boundary conditions, we find 
that they are transformed as follows: 


r)=0: f = g= 0 

T) = oo; f = C2 
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The last step is to transform the differential equations. Equations 
(3.11) and (3.12) become 


and 


pf + b prjf ' + g' = 0 


(3*29) 


pf 2 + | p-qff ' 


+ gf' = vf" 


(3.30) 


The above example illustrates the one-parameter group -theoretical method. 
Some further comments are now in order. 

In both examples, the steps are as follows: 

1. Define the group of transformation and substitute into the differential 
equations . 


2. Require that the differential equation be invariant. Relations among 
the constants in the transformation are obtained. 

3 . Eliminate the parameter of transformation to give absolute invariants, 
which will become similarity variables. 

4. Check the boundary conditions to see if they can be transformed into 
constant values, 

5 . Transform the differential equations. 

It is seen that up to step 4, no differentiation is needed. This makes 
the group-theoretic method an advantage over other methods. As long as in 
step 2 zeros for all the constants (e,g, , Oh., , . . ,a 5 in the linear group) are 
not obtained, the partial differential equations are transformable into 
ordinary differential equations. 

Step 4 in this particular problem needs not be put before step 5 J how- 
ever, in general, if the unknown function appears in the boundary conditions 
alone (e.g„, the diffusion equation with a boundary condition u(0,t) = U(t) 
and the form of U(t) is given), it ultimately may save effort to check the 
boundary conditions first as above. 

Equations (3«ll) and ( 3 . 12) can be reduced to one equation by solving 
Eq, ( 3 . 12 ) for g. However, we are illustrating the method and thus retain f 
and g in two separate equations. 
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Finally, comparison "between Eqs, (3,18) and (3=25) shows that Eq. (3.25) 
can he obtained, by putting a x = 0 in Eq.. (3.18). Thus, the spiral group is 
sometimes considered to be a special case of the linear group in which a x =0* 


3.2 MCJLTI PARAMETER GROUPS 

The method discussed in the previous article using a one -parameter group 
of transformation is usually applied to cases in which one independent var- 
iable is to be eliminated. If two variables are to be eliminated, the method 
should be applied twice, first reducing the variables by one and then apply- 
ing the technique to reduce the variables of the transformed differential 
equations again. The simplicity of the method characterizing this technique 
is then lost. Thus it would be valuable to have a method in which the number 
of variables can be reduced by more than one in a single application. Such a 
method based on multiparameter groups was presented by Manohar^ and later by 
Ames.^ We will discuss this method but limit the discussion to the method 
for reduction of the number of independent variables by two. Generalization 
to. more variables at a time can be made in a similar manner. Before proceed- 
ing with the nultiparameter group analysis, we shall first describe a special 
case where one-parameter groups may be employed to reduce the number of inde- 
pendent variables by more than one. 


3 . 2.1 One -Parameter .Group Method* 

Let a transformation group Fn be defined by 

a !_ a 2 _ o: 3 _ 

rn: x i = A x x , x 2 = A x 2 , x 3 = A x 3 

a i— 

Yj = A yj , (,j= 4 ,.. ,,n) 


(301) 


where A is the parameter. of transformation and the a’s are constants to be 
determined from the condition that the given partial differential equations 
be constant conformally invariant under this group of transformation. We con- 
sider a partial differential equation in which the independent variables are 
x 1; x 2 , and x 3 . 


In general, the group Pn can only reduce the variables by ope. However, 
in the special case in which a x = a 2 0, reduction of the variables by two 
is possible. Absolute invariants can be shown to be 


*Amesl considers the case in which a x =a 2 =Q as a separate case. However, since 
the form of absolute invariants in this case is the same as in r^4 below, 
this case is omitted. 
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( 3 - 32 ) 


(ax x +bx2 ) a 3/ Q: i 


and 


fj(n) 


(aXi+iDxa y i ^° 1 


where a and b are any arbitrary constants * 


(3.33) 


3=2»2 Two -Parameter Group Method 

Let r 2 i be a two-parameter transformation group 


A a-»_ 6*i_ ,ao 83— 

T 2 i: X! = A 1 x 1 , x 2 - B pi x 2 , x 3 = A 3 B P3 x 3 


a i Pi- 
J 3 - A 


( . . . ,n) . 


(3-34) 


Again, A and B are parameters of transformation and a's and p's are constants 
to be determined from the condition that the given partial differential equa- 
tions be conformally invariant under this group of transformation. Absolute 
invariants are found to be 


*n 


Xi « 3 /o;i X2 P3/Pi 


and 


(3.35) 




^5375^37^ 


Next, we consider the two -parameter group Tzz where 


(3-36) 


F22* Xj_ "* + Oi ^_A , 

a-vA. B-i- 
yj = e J B J y. 


3i— 

x 2 = B x 2 , 


a 3^P3— 

x 3 = e B x 3 


( j ^ } • • • } U. ) • 


( 3 - 37 ) 


Absolute invariants for this group are 


3 ^ 



( 3 . 58 ) 


■Q = 




«3 X 

“i %> PT 


P 3 


and 


f.i(n) = 




x x 

a i X2 Pi 


(3.39) 


A third group T 23 may he defined as 


Qh— o: 3 63B 

r23= *i = A x x , x 2 = x>> + PiB^ x 3 = A e x 3 


yj - A e 


(j=4, . . .,n) 


( 3 - 40 ) 


Absolute invariants for this group are: 


and 




Xc 


Xl 


0^3 
O'! 


.Pi 


Xs 


(3.41) 


fj(^) 



IA 



(3.42) 


Finally, let the group f 24 be defined as 


r 2 4* X x ~ X x t O^A 


, x 2 = x 2 + p x B, x 3 = e a3A e^ 3B x 3 


a.A PiB_ 
y . = e J e J y . 


( <3 -4 ) « « » , n ) 


(3.43) 


Absolute invariants for this group are: 




x 3 


oi 3 

,qs 


Xi 


M xs 


( 3 . 44 ) 
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and 


*J< 1 ) 


JX-, 


*1 


Pi 

:>31 


*2 


(3-^5) 


3.2.3 Example of a Similarity Analysis for a Three -Variable Problem 

Consider the unsteady, two-dimensional boundary-layer equations expressed 
in terms of a stream function 


chjr c) 2 ^ cty Su , SU S 3 \lr 

_j_ «. == 4- U — 4- y 1- 

dtdy 8y Sxdy 8x by 2 dt bx by 3 


(3-46) 


The "boundary conditions are 


y = 0: |t _ |t _ 0 

ay ox 

y = °°: ^ = u(x,t) . 


One -Parameter Group: 


Let I’ll Le defined as 


rn: t = A^f , x = A as x , y = A^y , 

4 = A^Sjr , U = A^U . 


(3.47) 


From the invariance of Eg. (3.46) we get 


1 1 

0:3 — g 0!i ^ 0!4 — OJ2 - Q-l ^ Q5 — CX2 “ 0!i 


(3-48) 


Reduction of the number of variables by two in a single step is possible 
only if a ± = a 2 , which gives 


a-. 


= 2 ai 


a 4 


a 


a > 


Old: — 0 


( 3 . 49 ) 
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Absolute invariants are: 


T) = - y - ■ (3.50) 

vat+bx 

and 

f(Tl) = ~ Z • (3.51) 

v (at+bx) 

Since U is a function of t and x, the only possibility for similarity to exist 
is to assume 


U = constant = U Q » 


(3.52) 


Using these variables, Eq* (3*46) can be easily transformed into an 
ordinary differential equation* 


3" t 


Is 

| T)f " + - ff" = 


(3.53) 


subject to the boundary conditions 


tj = 0: f = f 1 = 0 

TJ + 00 : f * = la 


Two -Parameter Group: 

Let Fax "be defined as 


r 2i : t = A a H, x = B Pl x, y = A as B p3 y 

t = A a4 B^ 4 t , U = A a5 B Ps U » 


( 3 * 54 ) 


From the invariance of Eq.. (3*46), we get 
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( 3 - 55 ) 


o% = -«i. 



a 4 = - 


1 

2 




P 4 - P 5 - Pxj p3 = 0 - 


Absolute invariants are : 






c 



X 


(3.56) 


Since U is a function of t and x, c cannot be a function of r|. Then the only- 
possibility is c to be a constant. Equation (3»46) is transformed into the 
following ordinary differential equation: 


ft + f« + | nf" - f' 2 + ff" + c(l-c) = 0 


( 3 - 57 ) 


subject to the boundary conditions 

T) = 0: f = f ' = 0 

T) = 00: f ' = 1 , 


It can be shown that for the other three groups mentioned earlier, all 
a’s and p's are zero. This means that similarity solutions do not exist for 
those groups. 
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CHAPTER 4 


GENERALIZED GROUP- THEORETIC ANALYSIS 


4.0 INTRODUCTION 

In the last chapter , the application of the simple group-theoretic method 
developed by Birkoff and Morgan-^ was treated. Extension of this technique 
to multiparameter transformation hy Manohar^ and Ames^ was also discussed. 

As it turns out, these methods are probably the simplest to apply of the 
principal techniques. Whereas other methods may involve solutions of differ- 
ential equations or other fairly complex mathematical manipulation, the group- 
theoretic method require straightforward algebraic procedures. One drawback 
of course, is that boundary conditions are not taken into account until the 
analysis is largely completed. However, checking for satisfaction of boundary 
conditions is a simple matter. 

_ O 

From the mathematical point of view, the method developed by Birkhoff 
and Morgan-^ considered two particular groups of transformation, namely, the 
linear and the spiral groups of transformation. Application of this method 
of similarity analysis to a given partial differential equation therefore 
answers only the question as to whether similarity solutions exist for these 
two groups. Since there is no proof that these two groups are the only two 
possible for similarity solutions to exist for a given partial differential 
equation it is still necessary to raise the question: Given a partial dif- 

ferential equation, what are all possible groups of transformation that make 
similarity solution possible? Are there groups other than the linear and 
spiral groups . 

To answer questions of this kind, we shall develop in this chapter a sys- 
tematic procedure in searching for all possible groups of transformation to 
a given partial differential equation. The procedure is based on Lie's theory 
of "infinitesimal continuous group of transformation" and his concept of con- 
tact transformation. Although both concepts were introduced by Lie in the 
latter part of the nineteenth century, there has been little application to 
the solution of nonlinear partial differential equations. The present tech- 
nique therefore can be considered as extension and application of Lie's the- 
ories of infinitesimal and contact transformations to the similarity solutions 
of boundary value problems. 

The technique itself, instead of the mathematical theories, is empha- 
sized in the development to follow. For a complete treatment of the theory, 
the reader is referred to the books by Lie9 - H and other references . 5 j 15 In 
the first section, definitions and theorems are presented briefly. This is 
followed by three examples illustrating techniques. The diffusion equation 
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is treated first. The goal of the analysis is to show how all the possible 
groups of transformation for a given partial differential equation (linear 
or nonlinear) may he found® By requiring that the given differential equa- 
tion he invariant under a "infinitesimal transformation/' the functional form 
of the so-called "characteristic function" can he determined. Once this is 
done* all possible groups of transformation can be obtained by the method 
developed in this chapter® 

The second example shows the application of this method to problems con- 
taining an arbitrary function that is to be determined as part of the analysis. 
The equation discussed is the steady* two-dimensional laminar boundary-layer 
equations. 

In the third example* the role of coordinate systems on similarity anal- 
yses is investigated by analyzing the Helmholtz equation in general orthogonal 
curvilinear coordinates* The problem encountered is the determination of the 
characteristic function mentioned earlier* By requiring that the Helmholtz 
equation be invariant under an infinitesimal transformation* a differential 
equation is obtained for the determination of the characteristic function. The 
equation contains the metric components (scaling factors) of a curvilinear 
coordinate system* The equation determines whether or not similarity solu- 
tions exist for a given coordinate system; or* given a group* what are the 
conditions under which the scale factors* h’s* should satisfy so that sim- 
ilarity solutions exist* 


4,1 DEFINITIONS AND THEOREMS 


4,1,1 DEFINITIONS 

A group is said to be continuous if* between any two operations of the 
group* a series of operations within the group can always be found of which 
‘the effect of any operation in the series differs from the effect of its pre- 
vious operation only infinitesimally* Thus* the transformation 


x* = a x x + a 2 y 
y* = b x x + b 2 y 


(4*1) 


is continuous if-af* a 2 * b 1; and b 2 are any real number* For example* if aio* 
a 20 * b 1Q * and b 20 are values of a a * a 2 * bx and b 2 * respectively* that carry 
(x 1 *y 1 ) into (x£*yj_)* a sequence of values of these parameters can be found to 
affect the same result* Each transformation can be made to differ from the 
previous one infinitesimally* 


4l 



The concept of infinitesimal transformations comes as a natural conse- 
quence of the definition of a continuous transformation group. An infinites- 
imal transformation is one whose effects differ infinitesimally from the iden 
tical transformation. Any transformation of a finite continuous transforma- 
tion group which contains the identical transformation can he obtained by 
infinite repetition of an infinitesimal transformation. 

Let the identical transformation be 


x x = ^(x,y,a 0 ) = x 
7 i = tU,y ; ac>) = y 


(4.2) 


where a 0 is a particular value of the general parameter a. Then the transfor 
mat ion 


xi = ^(x,y,a 0 +8e) 

yi = ■'K x ,y,a 0 +5e) , 


(4.3) 


where 5s is an infinitesimal quantity, defines an infinitesimal transforma- 
tion in a broad sense. A slightly more restricted definition will shortly be 
given based on the concept of 5s being an "infinitesimal" i.e., a quantity 
such that higher orders of 5s than 5s itself may be neglected in a given 
operation. 

Equation (4.3 ) can be expanded in Taylor series 'which gives 


and 


x i 


^( x ^yj.a 0 + s e) 


^(x,y,a 0 ) 


1 * V ^ a / a o ^ i \^a 2 / a 0 


yi 


i|r(x,y,a) = t( x ^y^ a o +5e ) 



Equations (4.4) and (4.5) can be written as 


(4.4) 


(4.5) 
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X + |Be + , . . 


Xi = 


C^-6) 


and. 


yi = y + iBe + ... (4.7) 

where 

% = i (, x >y> a o ) = fijf) (4.8) 

\ /Qjq 

and 

n = T|(x,y,ao) = ‘ (4.9) 

The expression /(x,y, & q ) and i|/-(x,y,a 0 ) in Eqs. (4,4) and (4*5) can he re- 
placed, respectively, hy x and y according to Eg,* (4,2). Now since a 0 is a 
fixed value of the parameter, a, both § and T] will 'be function of x and y only. 
Also, since Be is infinitesimal, we shall neglect higher order of Be and write; 


x x = x + |(x,y)8e (4,10) 

and 

7x = y + T](x,y)8e » (4,11) 


This transformation shall. "be defined as infinitesimal transformation. It can 
he shown that any one-parameter group, G x , contains only one unique infinites- 
imal transformation. The definition can he extended to a transformation con- 
taining any number of variables (hut only one parameter). 

Geometrically, an infinitesimal transformation transforms a point (x,y) 
to a neighboring point a distance 


n/ {(x 1 -x) 2 +(y 1 -y) 2 ) = nT (| 2 +T] 2 ) Be 


(4.12) 


in the direction 9 where 


cos 9 = £ — 


(4,13) 
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sin 0 


( 4 . 14 ) 


4,1.2 Representations of Infinitesimal Transformations 

Let f(x,y) "be a generally analytic function of x and y. The effect of an 
infinitesimal transformation 


xi = ^(x ; y,a 0 +5e) 


( 4 . 15 ) 


7i = t(x,y,a 0 +Se) 


( 4 . 16 ) 


on f(x,y) will he to produce the quantity f(x 1 ,y 1 ) which, upon expanding in 
Taylor series, 5 becomes 


f ( x i^yi) = f {^(x,y,a 0 +&e), y(x,y,a +6e)} 


= f { (x+|5e ) , (y+T]Se ) ) 


- f(x,y) + |f (I 3! + T) } 


, oe f 2 cr- 1 0 

+ or- (6 + 


+ W- ti n & * (5) i n -\ 


6x' n '“ 1 6y 


.. + ) 


( 4 . 17 ) 


Lie introduced the "symbol"* 


*It can readily be seen that 


Uf = (df (x-^y.J/da) 


and in particular, Ux = | and Uy = T| , 
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TT-P t j. „ 

uf = 1 Si + 11 d^ * 


(4.18) 


It should "be noted that this symbol itself is not a transformation.. It 
"represents" a transformation, however ; and is determined by it. 

It can be easily shown that the higher order term in the Taylor series 
expansion for f ( x x , y x ) can be written in terms of U as: 


l^f = (I 2 — + 2 1 Tj + ii 2 ^ } 

dx 2 dxdy By 2 


(4.19) 


{| n ( n )gn -i . <?t 

dx n 1 dx^dy 


+ + r,n 1 

+ ... + t) n ‘ 

ay 


(4.20 : 


where U 1 ^ represents operating on f by the operator 


,(x,y) + T)(x,y) jL 

dx dy 


(4.21) 


for n times, 


Equation (4,17) then becomes 


f Ui,yi) 


f(x,y) + yt Uf + |r4 U 2 f + .... 

X a c. % 


(4.22) 


By using this expansion, it will be possible to deduce the actual equa- 
tions of the transformation generated by a given infinitesimal transformation. 
Some examples will make this point clear. 


Example 4.1. Consider the infinitesimal transformation represented by 


-v + x 
J A T“ 

ax ay 


(4,23) 


x x = x + | (x,y)5e 


yi = y + Ti^y) 56 
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where 


and 


S(x>y) = -y 

T}(x,y) = X . 


We now seek expression for x x and Yi as one -parameter transformation 
groups. This can be done quite readily as follows: First take f = x, then 

using Eq. (4.23 ) , 


Ux = 

-y 

U 2 x = 

-X 

U 3 x = 

y 

% 

X 

11 

X 


Substituting into Eq. (4.22), we then get 


Se 5e 2 5e 3 5e 4 

Xi = x -ir y -2T x + 3T y + — x 


x{ 1 - §£-+ pt 


S€ 3 ^ S€ 5 , 

-y(8 e 


= x cos Se - y sin Se 


Similarly, if f = y, Eq. (4.22) then gives 


(4.24) 


(4.25) 
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5e 5e 2 5e 3 5e 4 

yx = y + x - 2I y - 51“ x + ij4~ y + * *■ * 


= x(5s 


3 1 5' 


+ y£l - §7- + |f- - .«« } 


= x sin Se + y cos 5e 


(4*26) 


The results correspond to a form of the rotation group: 


x x = x cos a ■*- y sin a 
yx = x sin a + y cos a 


for which a D = 0 in our analysis* 


Example 4*2 » Consider the infinitesimal transformation represented hy 


Uf 


If + <*7 


df 

Sy 


i * e a y 


and 


Xx = X + |(x,y)5e 


yx = y + n(x,y)5e 


where 

S (x,y) = c x x 

and 

T)(x,y) = c 2 y 

we proceed as before to get expressions for x x and y x : For f = x, we get 
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Uf 

= CjX 






iff 

= cf X 


(4,27) 




U 3 f 

* « W 

= cfx 



Therefore, Eg.* 

(4*22 ) 

gives 








Cn5e 

, c 2 5e 2 

cfSe 3 



Xl : 

= x + 

irr x 

4- X + 

d 0 

— S-j — x + « » .» 
7* 




= x(l 

+ 

-r ? 

X* 

. cf5e 2 
+ 21“ + 

,*) 




c< 

= xe 

jbe 

* 



(4*28) 

Similarly, 

for f 

= J, we get 







Jx 

cs5e 

= ye * 


(4*29) 


If we let ao = 0 in our analysis and define 

A = e a (Note, e a o + ^ e = e^ e ) 
then Eqs* (4*28) and (4.29) become 

x x = A Cl x , jx = A C2 y * (4.30) 

This is seen to he the linear group.® 

Example 4«_3» Consider the infinitesimal transformation represented by 

Uf - Cl | +W | . (4-31) 
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We may easily obtain : 


x x = x + c^Be 


( 4 . 32 ) 


and 


, C2&e , cfSe 2 

yi = 7 + y + — — y + * * ■ 


21 


= ye 


c 2 5e 


(4.33) 


This is seen to "be a case of the spiral group: 


x x = x + c x a 


yi 


ye 


c 2 a 




4.1*3 Functions Invariant Under a Given Group 

A function f(x,y) is said to he an invariant of a group if it is unaltered 
by the transformations of the group; therefore, for the function to be in- 
variant, we have 


fUi^yi) = f(x,y) 


(4*34) 


It follows that the function is necessarily invariant under an infinitesimal 
transformation* From Eq* (4*22), we have 


f (xi,yi) 


f (x,y) + Uf + 


Se" 

21 


U 2 f 


(4*35) 


But if f(x,y) is to be invariant, the following condition must be met: 


Uf = U^f = U 3 f = **» = 0 


(4*36) 


However, since the condition 


Uf = 0 


(4*37) 
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implies 


U 2 f = U 3 f = = 0, 


(4.38) 


Eg.* (4.37) is the sufficient condition that Eq* (4*34) is true for all values 
of x, y and Se« Hence we may state the following theorem* 


Theorem 4,1. The necessary and sufficient condition that f (x,y) he invariant 
under the group represented hy Uf is Uf = 0, 

In other words, the nec'essary and sufficient condition that f(x,y) he in- 
variant under a one-parameter group is that it he left Unaltered hy the in- 
finitesimal transformation of the group* 

To determine the invariant function, f, it is sufficient to solve the 
equation 



(M9) 


k 

The equation can he solved using the method of Lagrange in the theory of linear 
partial differential equations. The procedure in the case of three variables, 
x, y, and f, can he written as follows* 

To find the .general solution of 

p ! || + p * 3? - R (W,0) 


where P, , P 2 and R being functions of x, y, and f, solve 


dx dy df 

Pi P 2 " E 


( 4 . 41 ) 


If the general solution of this system is 


u i = c i and us = c 2 


( 4 . 42 ) 


*The system, Eq* (4*4l), has only two independent solutions. 
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then 


S^i(ui;u 2 ) = 0 


( 4 . 43 ) 


or 


Ui = <f>2 Ote) , 


(4.44) 


where <f>x and 9(2 are arbitrary functions of Ui, us, and us respectively, will 
he the general solution of Eq* (4»4l). 

As an illustration of this method, consider the equation 


xf + yf = xy . 
ox dy 

Based on the theorem, we consider the system 

dx _ d^; _ df 

xf yf xy * 

This system of equations has two independent solutions 


( 4 . 45 ) 


(4.46) 


xy - f 2 = ex 


and 


x _ 

— — cp ■» 

y 


Therefore, the general solution to Eq* (4.45) is of either of the following 
forms: 


^(xy-f 2 , | ) = 0 (4.47a) 

or 

~ = <(2 ( xy-f 2 ) (4,4Tb) 

u 

It is seen that application of this theorem to the solution of Eq, (4,39) 
for the invariant function, f, as a function of x and y forms a special case 
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as, R = 0. Equation (4,40) becomes 


df . bf _ n 

Pl S P2 ¥ ■ 


clearly, one solution of this equation is 


f = constant 


and we may write 


Ui(x,y,f) = f = constant 

The general solution of Eq, (4.48) is therefore (cf, Eq. (4,44)) 
U£ = <j> 3 (f) ~ 52(3 (constant) = constant 

where ua is a solution to 


dx _ dy 
Pi P 2 

We therefore get the following theorem: 
Theorem 4*2. To find the general solution of 


where % and rj being functions of x and y, solve 


dx _ dy 

t ~ n 

Let solution of this equation he expressed as 

fi(x,y) = constant. 


(4,48) 


(4.49) 


(4,50) 


(4,51) 


(*. 52 ) 


(*•53) 


(4.5*) 
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This function is the invariant function for the infinitesimal transformation* 
represented by 


Uf 


df 

cix 


+ T1 


df 

8 y 


V 


Since Eg. (4*53) has only one solution depending upon a simple 
stant, it follows that a one -parameter group in two variables, 
and only one independent invariant* 


Example 4.4* Consider the rotation group represented by 


where 


Uf 



+ x 


be 

8 y 


/ 


SU,y) 


-y and r)(x,y) = x . 


• Solution to the equation 


dx 

-y 


dy 

X 


gives 


x? + y 2 = constant* 


This is the invariant function of the rotation group* 


Example 4*5* Consider the linear group 


tt-p df , df 

Uf , CligJ+esrg- 


where | = c x x and T](x,y) = C2y* Solution to 


(4*55) 

arbitrary con- 
x, y, has one 


(4.56) 


(4*57) 


(4*58) 


*The function f is an invariant of the one-parameter group obtained from the 
infinitesimal transformation* 



then a necessary and sufficient condition is again 


Uf = 0 



6 j) 


■which, in the case of n variables, takes the form 


i 1 ( x i » * • • , x n ) Mr + **• + £n( x i>‘ * * > x n ) w - 0 


dx x 


bx 


•n 


(4.64) 


Following the same reasoning as in two-dimensional case, the invariant 
functions can be obtained by integrating 


dx x _ dxg _ dx 3 
li I2 I3 


— o * o — 


dX; 


■n 


(4.65) 


Since there are (n-l) independent solutions to Eq. ( 4 . 65 ) a one - parameter 
group in n variables has (n-l) independent invariants . The invariant func- 
tions are therefore 


fyn( x i, * * ° ^ x n) — c m> (m— 1, . » * ,n-l) 


(4,66) 


and are the solutions to the system of equations given by Eq. (4.65). 


4.1.5 Relationships Satisfied by Differential Equations Admitting a Given 
Group of Infinitesimal Transformations 

In this article a basic theorem on the determination of relationships 
satisfied by functions admitting a given group of infinitesimal transformation 
will be given. This theorem is comparable to the theorems of Morgan discussed 
in Chapter 3 . It is this theorem which lays the groundwork for the reduction 
of the number of variables in a partial differential equations. 

Consider now a function 


F 


/ 





d k yi 

d(x x ) k 



(4.67) 


the arguments of which, assumed p in number, contain derivatives of y^ up to 
order k. Such a function is known as a differential form of the k-th order in 
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m independent variables. Designate the arguments by z 1} 


> z p > e *§ 




Z 1 = X X 
z 2 = x 2 


z n- 




P-i 


y n 


S(xm_i) 


t. 


jjn 


Thus, Eq.. (4.67) can be written in a simpler form as 


(4.67a) 


F = F(z!,.,.,z p ) 


(4,68) 


The function F(zi,«*. ,Zp) is said to admit of a given group represented by 


Uf = li(zi,***,zp) —+... + ^(zx,.:.-^) 


6f 
8 z t 


(4.69) 


if it is invariant under this group of transformation. Therefore, the function, 
F, admits of a group if 


UF = 0 


(4.7Q) 


or. 



8 f 

+ | 2 SL- + 

oz 2 


+ i 


dF_ 
P Sz p 


0 . 


(4-71) 


It was shown in the preceding article that there are (p-l) functionally 
independent solutions, or invariants, to this equation. If they are denoted 
ty 


^m = fy&( z i.> * * * > z p ) = constant, m=l,...,p-l . (4.72) 
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Equation (4.71) must be satisfied, i.e., 


+ ... + h $L 

8z x . 8Zp 


= 0 . 


C^-75) 


Nov, if a change of variable is made in the function F, given by Eq» 
(4,68). from ( Zj_ , * • • , Zp ) to (ri x , . « *,ilp_x> Zp ) , ve then get 


F(z 1 ,», = ,Zp) - t(rji, • « » ,Tlp-i,Zp) 


(4.74) 


The condition given in Eq, (4,71) "will still have to be satisfied. Thus, 


up = U U + ... + lp^- 

dz x oz^ azp 

» M n 

- 51 V + 52 V + «»• + 5r) V = 0 , 

OZi 0 Z 2 ^ OZ„ 


(4,75) 


Since \|r is a function of t)x, « , • ,1p-i,Zp, the chain rule of differentiation may 
be used to get 


UF = | x + ... + _jl. Mp- i + djL. ) 

\o4i ozi orip-i dzx dz p az x y 

+ 6a (^lL^i + ... + 

Va^x a Z2 ^ z i 6z p ^ z 2/ 

“ 4 " O a a 

+ 1 /k_ hi + ... + 

\St]x Szp ari p _x a Zp a Zp 8z p / 


at 

Mi 


ix hi + 

V fcx 


+ | 


hi 


p a z 


o «> » 

L ••• + ip i n |t 

^1p-i \ Sz i ^Z p J &Zj 


= 0 


(4,76) 


57 



Since t)!, . . . ^p-i are invariant functions, all the "brackets equal to zero 
[cf. Eqs. (4.72) and (4.73)]> the following important conclusion then results: 


^t( T li,...,r]p-i,zp) 

Sz p 


( 4 . 77 ) 


This means \|r is independent of Zp. Equation (4.74) can then he written as 


F(z x , . . . ,z p ) = • • °>1p-i) = 0 * 


( 4 . 78 ) 


Therefore, if a differential equation 


F ( z x , . . . , Zp ) — 0 


is invariant under the infinitesimal transformation, it must he expressible in 
terms of the (p-l) functionally independent solutions of the partial differen- 
tial equations Uf = 0. 

Like Morgan's theorem given in Chapter 3> this result is the foundation 
upon which the technique given later in this chapter is based. 


4.1.6 The Extended Group Concept 

Consider the one-parameter group of transformation 


x x = sz$(x,y,a); y x = i|r(x,y,a) . 


( 4 . 79 ) 


Suppose y is regarded as a function of x. Then if the differential coefficent 
p(= dy/dx) he considered as a third variable, then under this group of trans- 
formations it will be transformed to p x where 


Pi 


dx x 


+ 5F 
M + M 

8x 8y 


P 


P 


X.(x,y,p,a) . 


(4.80) 


It can be easily shown that the general transformation 
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Xi = .^(-x,y,a); 


7i = a); Pi 


= x( x ;>y,P>a) (4*8l) 


form a group. This group is known as the extended group of the group given 

in Eg.. (4.79) * 

Now for a change 6e in the parameter a, we may express Eqs. (4.79) as 


x i 



S(x,y) 





" 1 " a » o 


(4.82) 


Yl 


x + 


5e 

il 


7( x ,y) 



8r) 

8x 



"f" oca 


(4.83) 


The transformed coefficient p x will he 


Pi = 


d y + if (M dx + ^ + - 

Se /dg \ 

dx + IT V 3x dx + ^ d y i + 


= p + 


§£(|l + f|l. |i\ p _ 

11 ^° x \oy d xj dy 


p + — £ (x^y^p) + ... 

11 


Therefore, the extended infinitesimal transformation given hy 


(4.84) 


x i 


is represented hy 


X + (5e)g, yi = y + (&e)r| 
Pi = P + ( & e)C 


UF 




8f 

8 P 


where 



(4.85) 


(4.86) 


(4.87) 
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Extension of this concept to higher order derivatives can he made hy the 
same reasoning. 


Example 1. For the rotation group 

Uf = 

it can easily he shown that the extended group of infinitesimal transformation 
is 


represented hy 


8f 

- y ^ + 


df 

ay 


(4.88) 


= X + |(x,y)5e, 

yi = y + Ti(x,y)Be, 

Pi = P + CU,y,p)Be (4.89) 


where 


i = -y 

T) = X 


and, from Eq.. (4.87), 


i 



8y dx/ 


P 


p 2 

ay 


1 + p 2 


(4.90) 


(4.91) 


The symbol of this extended group of transformation is therefore 


Uf = -y + x + (1+p 2 ) . (4.92) 

ax dy dp 

Example 2. For the linear group, the extended group of transformation can 
he represented hy 

Uf = cjx + c 2 y + (c 2 -c 1 )p |^ . (4-93) 

dx dv do 
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Example 3° For the spiral group, the extended group of transformation can he 
represented by 

Uf = c x ~ + c 2 y + c 2 p . (4.94) 

dx dy dp 

4.1.7 Contact Transformations 

The theories developed up to this point are not complete. If a given 
partial differential equation 


F — F(z 1 ,»..,Zp) - 0 (4=95) 

is invariant under the infinitesimal transformation represented 'by Uf, then 

UF = | x — + ... + £ = 0 (4.96) 

bz x y OZp 

and the differential equation can 'be expressed in terms of p-1 invariant 
functions. These invariant functions are solved hy using the system of equa- 
tions. 




(4.97) 


as discussed in article (4.l4). In case a particular group of transformations 
is given, as in examples in the previous sections, the function £ 1 ,... , £p in 
Eq. (4.97) are known. The invariant functions can then he found hy solving 
Eq. (4.97)« These functions are the similarity variables. However, the re- 
sults only give similarity solutions for this particular group. What is more 
important, is to derive the transformation groups and not specifying them. 

We shall now lay the groundwork for this type of analysis. 

In the present section. Lie’s theories of "contact transformations" are 
introduced. These theories make it possible to express the transformation 
functions, Si,«.»,£pj in terms of a single function, known as the "character- 
istic function." Basically, Lie’s theories of contact transformation deals 
with the transformation of a differential equation in a general, highly ab- 
stract, way. The abstractness and the complexity of the theories prevent any 
extensive discussion. Here, we merely state some of Lie’s theorems without 
proof. This is sufficient for purposes of application. 
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Lie 9 defined, a "lineal element" of the plane to he the ensemble of a 
point and a straight line passing through that point. In two-dimensional 
Cartesian coordinates, the coordinates of the lineal element consist of the 
coordinates, x and y, and the slope p.of the line. 

Thus any transformation in x, y, and p may he considered as a transforma- 
tion of the lineal element. 

Consider now a curve defined by y = y(x) having a slope dy/dx = p at the 
point (x,y). Let the variables x and y be transformed by 


Xi = X(x,y) 

7i = Y(x,y) . 


(tf- . 98 ) 


The curve is transformed into a curve having a slope p x at (x 1 ,y 1 ) defined by 


Pi 



dY + 3 y 

8 x By 


P 


}Sp 

ax ay 


x(x,y,p) . 


( 4 - 99 ) 


Equation (4-. 99) can be considered to be a transformation that along with Eq» 
( 4 , 98 ) transforms an element (x,y,p) into (x x ,y x ,p x )» Lie called this trans- 
formation an "extended transformation." 

A property of an extended transformation is that it transforms the dif- 
ferential equation 


dy - pdx = 0 


(4.100) 


into the differential equation 


dy x - p x dx x = 0 * 


(4.101) 


*The differential equation (4.100) is transformed into the differential equa- 
tion (4,101) under the assumption that the Jacobian of the transformation 
( 4 . 98 ) is nonvanishing. 
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We may say that every family of lineal elements (x,y,p) satisfying (4.100) is 
transformed by an extended transformation into a lineal element (x^y^Pi) 
satisfying (4.101). 

For three-dimensional space , x, y, z, an "element" is defined by five 
quantities, namely, x, y, z, p, and q where p and q, like y' in two-dimen- 
sional space, are the slopes of a straight line passing through that point. 

The family of surfaces containing a given element would be defined by a single 
relation between the five quantities, that is, by a partial differential equa- 
tion. 


Lie therefore defines a "contact transformation" on a transformation of 
the lineal elements of the plane which leaves the Pfaffiian equation 


dy - pdx = 0 (4.102) 

invariant; that is, a transformation 


*1 

= X(x,y,p ) 


7i 

= Y(x,y,p) 

(4.103) 

Pi 

= P(x,y,p) 



is a contact transformation if the Pfaffiian equation, Eq. (4.102), is trans- 
formed to 


hyi - PidXi = 0 


(4.104) 


after the transformation. 

O i 

The above definition of contact transformation is extended by Lie° as 
follows. 

When Z,X x , . » . ,X n ,Px, .. . ,P n are 2n+l independent functions of the 2n+l 
independent quantities z,Xj., . ... ,x n ,p 1 , . . ». ,p n such that the relation 


dZ - Pj_dXj_ = p (dz-pqdxi ) * (4.105) 


*Pj_dXi PqdXq, i.e« ? the repetition of indices will indicate summation over 

all values of the index. 
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(where p does not vanish) is identically satisfied, then the transformation 
defined hy 

Z = 

Xn = X n (z,x 1 , . . . ,x n ,Pi, . . . ,p n ) 

?n = Pn( z > x i j * “ • ; x n.>P:u • • • jPn) (4.106) 

is called a contact transformation. 

The contact transformation, Eq. (4.106) will transform a partial differ- 
ential equation in z,x 1 , . . . ,x n ,p 1 , . « . ,p a into one in Z,X x , . . . ,X n ,P XJ . . . ,P n 
and also the solution of the first partial differential equation into the solu- 
tion of second. 

In the next article, the property of a contact transformation, Eq. (4.105) 
will he used to derive expressions for the transformation functions in terms 
of characteristic functions. 

4.1,8 Contact Transformation in Terms of the Characteristic Function 
First, let us consider the infinitesimal transformation 

Z = z+5ef;, Xq = xq+-8e|q, Pq = pq+Sexq (4.107) 

where £, and jtq are functions of z,x x , . , . ,x n ,p x , . . . ,p n and e is a quantity 
so small that its square and higher powers may he neglected. We now wish to 
find the functions £, |q, jtq in terms of a single function of z, xq, pq, known 
as the characteristic function. 

From the definition of a contact transformation 

dZ - PqdXq = p(dz-p i dx i ) (4.108) 

we get 


64 



(4.109) 


az , Sz , dz 3 „ /dX-s , BXi „ bx± 

d2 * 51 dXi + 51 dPl ■ Mar dz + S£ ^ + 5 ; dp r 


Sz 


= p(dz-p i dx i ) 


Equating coefficients' of dz, dxp, and dp^ on loth sides, we get 


az 
a z 

az 

Sp r 

az 

ax,. 


p 

1 


p 4 


axj 

1 ax-n 


= 0 


-PPv 


(r=l, o « . /n) 


(4.110) 


For the infinitesimal transformation defined in Eq., (4.107), Eqs. (4.110) become 


K 

az 


at 

i 


a 




SL..p. ik 
a Xr 1 a Xr 


0 


-*r = -dp r 


(4.111) 


where the relation p = 1+ecr is used (a, arbitrary). Now, p.j. is independent of 
x r , and z. Therefore, if we define a characteristic function W such that 


W 


Pi^i 


- i . 


(4.112) 


Equations (4,111) will yield 
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a 


6w 

6z 


6W 

^P r 

6x r 


5 ir^i = 6r * 

-a Pr + n r = p r + rt r . 


We then get 


and 


i 




- w 


3w dw 

a^ r a z 


a = 


dw 

bz 


(4.113) 


(4.114) 


(4.115) 


Nextj consider the infinitesimal transformation 


Zi 

= Zi + Bemi(z v ,x^,pJ) 


x k 

= X3£ + Bea k (z v ,xj 1 ,pJ) 


Pk 

= Pk+Se4( z v>V p i 

(4.116) 


■where m-^ Cfc and are functions of x x , , . . ,x m ; z ± , . . . ,z n ; and , . • • ,p^ 

(= dz n /dx m ) and e is a quantity so small that its square and higher powers 

may he neglected. We now want to express the functions m±, and 

(i=l, • • « ,n;k=l, . „ „ ,m) in terms of ¥ 1 (i=l ; .., ; n), the characteristic functions. 

From the definition of contact transformation we have 

p k dX k = Pi v ( dz v“Pj dx |j.) (4.117) 

the second equation in Eq. (4.111) the term 


dZ-i 


*By adding to both sides of 
^V± 


dpi 


^i' 
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we get 



SZi 

+ ^ dx x 


+ 


bZ ± 







+ 




= P iv (^v-P^) » 


(4.118) 


Equating, coefficients of dz V} dx^, 


and dp^ we get 


. pi 

8z v k 8z v 


^i.pi^ 
k Bxp, 


8Zi 


,i ^ 




iv 


" p iv p ^ 


0 . 


(4.119) 


For the infinitesimal transformation given in Eq.„ (4„ll6), Eqs. (4.119) hecome 


(“fpk 0 ^) 

OZy 

= %v 


A- (mj-pio^) 

= “Pp|Piv + 


(“i-P^k) 

6Pp 

= °y,6iy 

(4.120) 

where the relation = 6^ v +eRj_ v is 

function Wj_ as 

used. Let us now define the 

characteristic 

Wi = 

Pk a k “ m i • 

(4.121) 
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Equations (4.120) then hecome 



-p^iv + *s 


a ia 5 iv • 


(4.122; 


From these equations ve then get 


— pj - Wi 


hpv.^ 

8z v P 8x,, 


(4.123) 


Having expressed mp, Q^, and in terms of the characteristic functions Wp, 
we now prove a very important property of mp and 0 ^. 

For n=l, Eqs. (4.123) reduce to Eqs. (4.114) which are expressed in terms 
of a single characteristic function. 


For n > 1, however, there is a very important property involving m^ and 
0^.45 Now, differentiating the third equation in Eq. (4.120) with respect to 
pP, we get 


JL JL 


( ®i" p k Q k ) 


s P g < ' &Iv0 *' ) 


ijb ^ § • 


(4.124) 


The third equation in Eq. (4.120) can he written as 


- (mi-pi^) 


-SipCHcr 


( 4 . 125 : 
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Differentiating Eg.. (4.125) with respect to p v , we get 


dpnSpa 


(mi-pja k ) = -5 i p 


da c 

$vl 


(4.126) 


Since the left side of Eqs. (4.124) and (4.126) are the same, we can 
equate their right sides and get 


iv apg 

If we set i=v^p, Eq. (4.127) then gives 


Hi = 5> da a 


!P TV 




(4.127) 


^ = 0 . 

Spg 

This means Op is independent of p§. 

Next, the third equation in Eq. (4.120) gives 


Smi 


SpJ ~ 

a k — - - S 
k SpV 


iv^p 


(4.128) 


- Q«k®kp ” a p ~ °p. “ 


a,. = 0 


(4.129) 


This equation means mi is independent of p% 

As a result, the infinitesimal transformation defined in Eq. (4, 116 ) 
should he modified to : 

Zi = zi + 5emi(z v ,x ia ) 

Xk = % + 5 «%( z v; x |i) 

Pj = Ej. + B e «l (z v ,^,pV) (4.150) 

where mi, a k , and n k are defined in Eqs. (4.125). Equations (4.125) also sug- 

gest that the dependence of Wi on p^ is linear— a very far-reaching and im- 
portant conclusion which will he made use of later. 
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4,1.9 Contact Transformation of Higher Orders 

Let us now consider the extended infinitesimal contact transformation 
defined by 


Zi = z 1 +.-8em 1 (z v ,x^,pJ) 

X k = x k + -6ea k (z v ^,pj) 

P k = p i + - 6 ert k (z v>VPj) 

P jk = p jk + &ert jk*‘ z v^ x |i-’P(i^(Il!s) 

= p jU + ;B6Jt ^ (z v^n^ p ^ p ^s^ p ^t)' (4 ' 1 5 1 ) 


where a k and mj_ are independent of p 


v 

d 


for i > 1, 


and 



dz-j_ 

dxjdx k 


etc. 


(4.132) 


We now want to express it 


W, 


The functions mj_, a k 


jk 


in Eqs. (4.123). 


and it 
and 


i 

«Jki 


as functions of the characteristic function 
it k are expressed in terms of W-j. and given 


By definition, 


dB 


k 


" P jk^j 


( 4 . 133 ) 


Substituting the infinitesimal contact transformation from Eg.. (4.131) into 
Eq. (4.133), we get 


d 4 = p jk da j + 4k dx j • ^* 1 ^) 

Subtracting both sides of Eq, (4.134) by the quantity Ctjp|j k dx s , we then get 

d (k k “Pjk Q: j ^ = 4k dx <] " a j p sjk dx s * (4.135) 
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The second term on the right side of Eg. (4.135 ) should he dropped since in 
derivatives of order higher than the second are omitted. We then get 






or, 


in terms of aj and it 


i 

k" 



+ pS 

a Xj a Zs J 


+ 


M -rvV 






(4.136) 


(4.137) 


To express it^ i n "terms of the characteristic functions Wj_ , we have to put 
expressions for it^ and from Egs. (4.123) into Eg. (4.137). 

Similarly, for the third-order function we get 



d 

dxi 



1 a ) 
Jkt V 


(4.138) 


or, 



(4.139) 


In order to express in terms of the characteristic functions, we have to 

express and a % in terms of W^, as in Egs. (4.137) and (4.123), and sub- 

stitute into Eg. (4.139)* 

The same principle can he extended to higher order, e.g. , 


jkim 




Ujki-Pjkff 0 ^ ) 


(4.140) 
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4.2 SIMILARITY ANALYSIS OF DIFFUSION EQUATION 

We shall now look at the diffusion equation from another point of view, 
namely, the searching for all possible groups of transformation that will re- 
duce the diffusion equation to an ordinary differential equation. In applying 
such a technique to a given differential equation, it may turn out that for 
some or all of the groups other than the linear and spiral groups, the boundary 
condition cannot be transformed although the partial differential equation can 
be transformed into an ordinary differential equation. For such cases, we are 
at least assured that the groups of transformations that remain are the groups 
possible for the given boundary value problems. A similarity analysis of the 
diffusion equation from this point of view is apparently not covered in the 
literature. The one -dimensional form of the diffusion equation in rectangular 
coordinate is chosen because of its simplicity. Extension of analyses to equa- 
tions expressed in other coordinates can readily be made. 


4.2.1 Analysis 

Consider the diffusion equation 


du _ 3 2 u 

St ciy 2 


0 


(4.142) 


on which an infinitesimal transformation is to be made on the dependent and 
independent variables and derivatives of the dependent variable with respect 
to the independent variable. The infinitesimal transformation is 
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t> = t + 8eg(t,y,u,p,q) 

y’ = y + B eTl(t,y,u,p,q) 

u' = u + 8e£;(t,y,u,p,q) 

P* = P + BertiCt^y^p^q.) 

q. 1 = q. + &eit 2 (t,y,u,p,q) 

P^2 = P22 + Sejr 22 (t,y,u,p,q,p 11 ,p 12 ,P22) ( 4 . 143 ) 


where , in terms of the characteristic function, W, 


dW 

a P 


aw 

Sq. 


i = 


-*1 = 


op Scl 

dW , bW 
bt V bu 


-Jt 2 


bW . „ bW 
’ & 


-^22 = 


+ 2, |fW_ + 4 , 




5F 1 4 s? P “ ssf. 


>w + 4 fv_) + p?2 afw + 2pi2P22 afw 

op opoq. 


+ 2p 22 f q^*- 

V^q ^q 


2 a w ^ „ aw 


. ^ VV , „ 

+ P 22 — + P 22 

aq 2 au 


(4.144) 


The characteristic function, W, is a function of t, y, u, p, and q» We note 
that 


au au a 2 u 

p " at' a a y ' P22 ~ ay2 


, P12 = — . (4.145) 

at ay 
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It is shown in article 4.1.3 that the necessary and sufficient condition 
that a partial differential equation F ( t , y .u.p , q . Pjo .P pp ) = 0 invariant under 
the group of transformation represented by Uf is UF = 0 which for the diffusion 
equation, is 


U(p-vp22 ) = 0 


(4.146) 


or, expanding the expression hy employing the operator U: 






3u” 




“5T 


+ rtll iU + * 22 iU 


Ml — - — — t Jt22 — — — + rti2 — — - = 0 

dpil 3p 22 Spi2 


(4.147) 


where the parenthesis represents the differential equation 


P - vp 2 2 


(4.148) 


Carrying out the operation in Eq. (4.147) yields: 


rti - VJt22 = 0. 


(4.149) 


Upon substituting expressions from Eq. (4.144) into Eq. (4.149) yields 


8W 8W . 3 2 W . 0 8 2 W . 2 

- — - v — + v — + 2vq — — + vq^ 

at ou Sy 2 8y8u 


6 2 W 

8u 2 


o„ 8 2 w „ A 0 „ „ 8 2 w 

+ ^-Pl2 v \ V + ^Pi2q yy - 
oyop auop 


v + 2p 


6 2 W 

dydq 


+ pqft- 

ouoq 


+ vpf 2 


8p £ 


2 8 2 w _ 8 2 w p 2 8 2 ¥ 8w 

~ + 2pi2P 5ST7a7 +p 7 " 0 • 


(4.150) 


Equation (4.150) is seen to be a linear partial differential equation in 
W(t,y,u,p,q). 


Since W is not a function of P 12 , the coefficients of the terms involving 
P 12 and p 2 2 should be zero. We then get 
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0 


(4.151) 


6 2 W 

8p 2 


6 2 ¥ 6 g ¥ p 6 2 ¥ 

6y8p ^ 8udp v cipdq 

Equation (4,151) indicates that W is a linear function of p. 
write 


(4.152) 

Thus we can 


¥ = Wx(t,y,u,q.) + p¥ 2 (t,y,u,q) . 


(4.155) 


Substituting this form of ¥ into Eq. (4.152 ), we get 


c3Wp 8W 2 p 8W 2 

aT + + v aT 


0 . 


(4.154) 


Since ¥ 2 is not a function of p, the coefficient of p in Eq.. (4.154) must be 
zero. ¥e therefore obtain two equations, namely, 



0 


(4.155) 


Ik + ^ = o . (4.156) 

8y ou 


Equation (4.155) indicates that ¥ 2 is not a function of q, i.e,, . 

¥ 2 = ¥ 2 (t,y,u), and so Eq. (4.156) can be broken into the two equations, 


and 


1^ = 0 (4.157a) 

ay 


^- = 0. (4.157b) 

du 


This means ¥ 2 is independent of both y and u and, as a result. 


¥ 2 = ¥ 2 (t) 


(4.158) 


and the characteristic function now takes the form 
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w 


Wi(t,y,u,q) + pWa(t) . 


( 4 , 159 ) 


Patting this form of W into Eq. (4.150), we get 


5W-> „ 6W 2 , 

— — =4. - p ==» + 

at at 


V 


afw 

dy 2 


+ 2vq_ 


aggi + 
ayau 


vq 2 


a 2 w 1 

au 2 


+ 2p 


afw^ 

dydq. 


+ pq 


^ 2 Wi + 

a u aq 


pf. a 2 w! 

v aq 2 


= 0 . 


(4.160) 


Since both W x and. W 2 are independent of p, Eq. (4.160) can be separated 
into three equations, corresponding to the coefficients of p°, pi, and p 2 . 

We then get 


aw T 

v o 2 W x 

+ 2vq --- 1 - 4 

• vq 2 

= 0 

( 4 . 161 ) 

at 

ay 2 

dydu 

au 2 


pi; 

aw 2 

+ 2 + 1 

a 2 w x 

= 0 

( 4 . 162 ) 

Jr 

at 

a y aq 

a u aq 



p 2 : 

a 2 w x _ 0 

0 


( 4 . 165 ) 


aq 2 


From Eq. ( 4 . 163 ), W x is linearly dependent on q, therefore, it becomes 

W x = W u (t,y,u) + W l2 (t,y,u)q . (4.164) 


Putting this form of W x into Eq, (4„l62), we get 

_ + 2 ^ + q = 0, (4.165) 

dt 6y au 

Both W 2 and W l2 are independent of q, therefore, Eq. (4,l65) becomes 

- S™ 2 + 2 ^ W i g = 0 (4.166) 

dt Sy 

6w l2 _ 
au 
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( 4 . 167 ) 



From Eq. (4.167), ¥ 12 is independent of u, Also, since ¥ 2 is a function of 
t only, Eq.. (4.l66) shows that ¥ l2 depends linearly on y, i.e„, 


¥ l2 = ¥ l2l (t) + ¥ l22 (t)y 

Equation (4.162) then becomes 


(4.168) 


- + 2 ¥' 
dt 


122 


= 0 


(4.169) 


¥e will make use of this equation later. 

The characteristic function, ¥, now becomes 

¥ = ¥ 11 (t,y,u) + {¥ l2 i(t)+¥ l22 (t)y}q + W 2 (t)p 

Putting Eq. (4.164) into Eq. (4.l6l), we get 


(4.170) 


_ dWn. f q + V 

dt \ dt dt J 


+ 2vq + vq 2 = 0 

' ' ciu 2 


Zy 2 ‘ 8y8u du‘ (4.171) 

Since ¥ X1 , ¥ 12l and ¥ 122 are independent of q, terms with different powers of 


q. are grouped and their 

coefficients are put equal to zero* 

Three equations 

are obtained: 

4 : -tf + v ch! 0 

(4.172) 

qi: 

_ jjjiai _ Mfgg y + 2v - g - 0 

dt dt dydu 

(4.173) 


p 8 ¥i 2 r\ 

q. • ■ = o » 

(4.174) 

Equation (4.174) shows 

that ¥ X1 is linearly dependent on u. 

Therefore, 


Wn = ¥xn(t,y) + W ll2 (t,y)u . 

(4.175) 
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Equation (4,173) then gives 


_ dW i2x _ dW iaa y + 2v = o . 

dt dt 8y 


(4.176) 


Therefore, ¥ 11 2 can "be written as 


¥ x 12 = Wxx2x("t) + Wx 1 22('t)y + W 1 i23(t)y 2 


Equation (4.176), then becomes 


"dft^ ~ y + 2 vW 112 2 + 4vyW xl 23 = 0 


Since all the ¥’s in Eq. (4,178) are independent of y, we get 


d ^.ia 1 4 . 2vV 1 x 22 

dt 


= 0 


+ 4v¥ 112 3 
dt 


= 0 


Putting ¥ 1X into Eq. (4,172), we get 


tylJUL + v 6 2 ¥xxi . 


dt 


Syg 


+ dW naa y + ffl ii23 yS. 
dt dt dt 


u 


+ 2v¥ 1123 u = 0 


Equation (4,l8l) can he separated into: 


u° : - 


d w m + v 6 2 ¥ 11X _ 0 


St 




(4.177) 


(4.178) 


(4.179) 

(4.180) 


(4.181) 


(4,182) 
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0 


(4.183) 


u i ; _ aw 112x + 2 vw 1123 = 

dt 

dW lx22 0 
dt 

d Wjj_g_3 _ o 

dt 

From Eq.s, (4„l84) and (4.185), 

Wn22 = Cl 

and 

^1123 = c 2 * 

From Eg., (4 ,183), 

w ii2i = 2vc 2 t + e 3 * 

From Egs. (4,179) and (4,l80), 

¥ 121 = 2vc 1 t + c 4 

^122 = 4vC 2 t + C 5 a 

From Eq., (4,l69)> 

W 2 = 4vc 2 t 2 + 2c 5 t + c 6 „ 

The final form of the characteristic function is therefore 


(4.184) 

(4.185) 

(4, 186 a) 
(4.186b) 

(4.187) 

(4.188) 

(4.189) 

(4.190) 
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(4.191) 


where W llx (t,y) Is any function satisfying Eq. (4.182), i.e., 


d w ni _ v . ^ 1.11. _ o 
at v ^y2 “ 0 


(4.192) 


The characteristic function, W, given in Eq. (4.191), will now be used 
to determine the absolute invariants. 

4.2.2 Absolute Invariants and The Transformed Equations 

With the characteristic function, W, obtained as in Eq. (4.191), we now 
make use of the general theory to find the absolute invariants. From Eq. 
(4.65), the following relations are obtained: 


dt _ dy _ du 

I 4 “ 5 (4.193) 


where the transformation functions |, r), and £ can be obtained by putting into 
Eq. (4.144) the characteristic function, W, given by Eq. (4.191,,). Equa- 
tion (4.193) then becomes 


dt _ dy_ 

4vc2t H +2c 5 t+c 6 2vc 1 t+c 4 +(4vc2t+c 5 )y 

= du . (4.194) 

-Win (1 ,J ) - ( 2 vc 2 t+c 3 +c 1 y+c 2 y 2 )u 

The number of possible groups are large, due to the fact that all c's 
are arbitrary and W xxi (t,y) is an arbitrary function satisfying Eq. (4.182). 
Therefore, we investigate a few special cases of the parameters. Other 
groups can be obtained in a similar manner. 
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Case 1. W ni (t ; y) = c x = c 2 = c 4 = c 6 = 0 


Equation (4. 19*0 becomes 


dt _ dy _ du 
2c 5 t c 5 y -c 3 u 


(4.195) 


The two independent solutions to Eg., (4.195) are 


and 


y_ 




constant 


~ constant: a — * 

t a 2 c 5 

According to the theories in the preceding articles, the diffusion equation 
can he expressed in terms of these two invariants, i.e., 


and 



f (t) ) 



(4,196a) 


(4.196b) 


The diffusion equation is then transformed into an ordinary differential equa- 
tion 


vf" = at - - rif 1 . 

2 

The transformation are seen to be the linear group of transformations. 


Case 2. W li;L (t,y) = c x = c 2 = c 4 = c 5 = 0 
Equation (4.194) then becomes 


dt du 

c 6 c 3 u 


(4.197) 
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and 


y = constant. (4.198) 

Following the same arguments as in Case 1, we get the absolute invariants as 

n = y and f (t) ) = ; (p = - Sa ) 

e Pt c 6 

&nd the diffusion equation is transformed to 

vf" - pf = 0 . 

The transformations are seen to 'be the spiral group. 

Case 5. W i;L1 (t,y) = c 2 = c 3 = c 4 = c 5 = c 6 = 0 
Equation (4,194) becomes 


.. jy = 

2vc x t 


du 


-ciyu 


and 


(4.199) 


t = constant* 


(4.200) 


The absolute invariants are 


tj = t and f (tj ) = 


u 




and the transformed equation is 


2r\f • + f = 0 . 
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Case 4, W i;L1 (t,y) = c x = = c 5 = 0 

Equation (4,194) "becomes 


dt _ dy du 

C 6 C4 -C3U 


The absolute invariants are 


T] = y - t and f 
c 6 



(4,201) 


%- 

The diffusion equation becomes 


vf" + ^ f 1 + ^ f 
c 6 C 6 


0 „ 


The above four cases are examples of cases where the solution of Eq, 
(4,194) is straightforward, i,e,, the two independent solutions can be solved 
by simple pairing of equations. The following cases are those in which the 
two solutions have to be solved in a sequence of steps. 


Case 9 , W xll (t,y) = c x = c 3 = c 4 = c 5 = c 6 = 0 
Equation (4,194) now becomes 


dt _ dy 
4vt a 4vty 


du 

-(2vt+y 2 )u 


(4,202) 


The first of the two equations gives the solution 


y = 

t 


1 _ 


(4,203) 


Next, replacing t in the last two terms of Eq, (4,202) by k x y [based on Eq, 
(4,203)], we get 


(2vk 1 yty 2 ) dy _ du 
4vk 1 y 2 u 


(4,204) 
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The solution to Eq. (4.204) is 


uy 1 / 2 e y/4vk x _ cons ^ an ^ 


or, using Eq. (4.203) again, 


uyl/2 


y 2 /4vt 


= k 2 


(4.204a) 


According to the applicable theorems, the absolute invariants are, from Eqs, 
(4.203) and (4.204a) 


r| = 21 and f ( T) ) = uy 1 / 2 t 

U 


The diffusion equation is then transformed into 


T] 2 f" - rjf ' + | f = 0 . 


Case 6. W xll (t,y) = c 2 = c 3 = c 4 = c 5 = 0 
Equation (4.194) becomes 


dt _ dy 
c 6 2vc x t 


du 

-c x yu * 


The first two terms give 


y - 



= 


(4.205) 


where k x is the constant of integration. Combining the first and the third 
term and making use of Eq. (4.205), we get 


(kit + & t 3 ) 

ue C ® 


= k 2 


(4.206) 


where k 2 is the constant of integration. Using Eq. (4.205), we get 
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(yt - fgi t 3 ) 


ue 


=, k 2 


(4.207) 


Equations (4.205) and (4.207) give the invariants. 


n = y - ^ t 2 
c 6 

and 


f(il) = ue 


^ (yt - t 3 ) 


c 6 


The diffusion equation is transformed to 


f" + rif = 0 . 

c 6 v 

Case Wiii(t ,y ) = c 2 = c 3 = c 4 = c. s = 0 

Equation (4,194) becomes 

dt _ dy _ du 

2c 5 t 2vc 1 t+c 5 y -Ciyu ° (4.208) 


By following the same steps as in the two previous cases, the invariants are 
found to he 


and 



- iv£LL & 

C 5 


f(i) 




The diffusion equation is then transformed to 


vf" 



0 . 
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In all the above cases, W ui (t,y) was taken to be zero. This is not 
necessary as we shall show in the following two cases. 


Case 8. W llx (t,y) = va x t + - a^ 


C x = C2 = C 3 = C 4 = Cg — 0 * 


It can be shown easily that this functional form of W lxl (t,y) satisfies 
Eg.. (4.192). For this case, Eq.* (4.194) becomes 


dt _ dy _ du 

2C = t ‘ C 5 y ‘ _ vait . | ai y2 ' 


Using the same method as in previous cases , we find 


and 



f(n) = u + a hf vt *y 2 > . 

4c 5 


The diffusion equation is transformed to 


Vf" + | T|f ’ = 0 


1 2 

Case ,9. W lxl (t,y) = a 0 +va x t + - a x y 

2 


c x = C2 = C3 — C4 - Cg - 0 . 

The only change made in this case is the addition of a constant term, a Q , 
to W xxx (t,y). The invariants in this case are 


T) 
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and 


f ( t) ) = u + ^(.f.x t+ y s ), 

4c 5 


a Q 

2c c 


in t 


The diffusion equation becomes 


vf" + - T]f ' - ^2_ = 0 . 

2 1 2c 5 


It is seen that as a result of the additional constant term, a 0 , one more term 
is added to f(r]) and the transformed equation, as compared with Cabe 8. 


4„3 SIMILARITY ANALYSIS OF STEADY, TWO-DIMENSIONAL, LAMINAR, BOUNDARY-LAYER 

EQUATIONS 

This article is the second application of the theories given in article 
4.1. The example treated in the previous article concerns the diffusion equa- 
tion. The present article differs from that example in that an unknown func- 
tion is involved in the differential equation which has to he determined. This 
factor introduces complexity in the process of searching for the possible groups 
of transformation. The method developed in this article can he used for simi- 
lar problems. 

We shall see that the houndary-layer equation may he transformed from a 
nonlinear partial differential equation to nonlinear ordinary differential 
equation by groups of transformation other than the linear and spiral groups. 
However, in transforming usual boundary conditions, the linear and the spiral 
groups are the only two groups possible. 


4.3.1 Infinitesimal Transformation and the Characteristic Function 

The laminar, incompressible boundary- layer equation, in terms of the 
stream function i]/, can be written as 


cty d 2 t _ d\|/_ d 2 i)r _ u d- u e + d 3 -^ 
dy dxdy dx dy 2 e dx dy 3 


With the boundary conditions 

y = 0: 


cfy 

dx 

dt 

8 y 


^ = o 
°y 

u e (x) 


(4.209) 
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where all the quantities are in dimensionless form "by the transformation 


x = 


L 


T 


J-l 


* = 


E_ N^Re 


LU r 



Re 


Mi 

V 


where x', y', ijr r and Ug refer to dimensional quantities and L and U D are 
reference length and velocity. 


Equation (4.209) can he put in a shorthand form as 


P2P 12 - P1P22 = i + P222 


»- (4.210) 


where 


Pi 


a? 


wS Pl2 

dy 


8x8y 


, etc . } 


and 


= Up 


du e 

dx 


(4. 21i; 


" Wow, let us denote the following differential expression hy F: 


F - P 222 + i - P 2 P 12 + P 1 P 22 » (4.212) 

The equation F = 0 will he invariant under the infinitesimal transformation 
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x' = x*-6ea!x(x,yy|rjPijPa) 
y* = y+.-aeCfe(x,y,T|r,Px,B 2 ) 
t’ = if +-5em(x,y,ijr,Pi,P2) 

Pi = Pi + 8 ejt!(x,y,i|r,Pi,P 2 ) 

Pi = P2+8ejta(x,y,t,Pi,P2) 

P 12 = Pi 2 +- 8 e«i 2 (x,y,i|r,p 1 ,p 2 ,Pii,Pi 2 ^P 22 ) 

pi 2 = P 22 +-Be« 22 (x,y,t,P 1 ,P 2 ,p 1 i,Pi 2 ,P 22 ) 

v- 

pi22 = P222 +'8ejr222(x,y,l|r,Pi,P2 jPii>Pl2jP22>Plll> • • • >P222 ) (4.213) 

UF = 0 (4.214) 



The next step is to express all the transformation functions, ct 1} (X 2 , etc., 
in terms of the characteristic function, ¥. The functional form of ¥ is then 
determined "by Eg.. (4.216). Now, from Eqs. (4.114), (4.137) and (4.139)* 


cc 2 

*1 

Jt2 

m 


dpi 

d¥ 


Sp 2 


'2 

8w 

8w 

SF Pl 

" s 


8 ¥ 

p 2 

8i|r 

“ ^y 



(4.217) 
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and 


-^12 


JfW 

dxdy 


+ P2 


Sx^Jr 



+ p 2 



J. -n j. „ 

+ Pill — + P2 

\0Pidy 



+ 



8w 

cty 


P2 


8 S W \ 

^P2^ly 


+ 



d 2 W „ J. „ 4. 

Pi + — — Pll + 

dtdpi dpf Sp 2 Sp! 



+ 



+ 


c^W _ L 8 2 W _ 8 2 W 

r~r — Pi + 7 — r — Pii + 7-r 
oi|fop2 opiop 2 dpi 



-Jt22 


dfw 

dy 2 


+ 2p 2 


8 2 W 

dy^|r 


+ PI 


8 2 w 

8\|r 2 


+ 2p l2 


8 2 w 

^y^pi 


+ p 2 


8 2 w 

<tydpi. 


2pssl 


^ydp 2 



S P ? 


2P12P22 


Spldp 2 


+ 


2 

P22 


8 2 w 

dpi 


+ p 22 


8w 

cfy 


( 4 . 218 ) 


( 4 . 219 ) 
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-„ 222 = . ^ . as* p 2 . ^ _ Is* 

0X2 OZ op^ p dp bc 

23L + 2 Pa J2L 

Sx§ Sxi^2 


3Pbo 


+ pi 


J2»_ + 2 Pl2 fJ2L_ + pg . . . ^ 

SX2^Z 2 VSxfSp! Bx 2 3zBp^ 

2p 2 2 ( - i 3 -W... + p 2 + P f 2 J!w__ + 2p l2 p 22 ^!w 

lSx 2 S P 2 Sx 2 c)zS P 2 y Sx 2 S P f Sx 2 ^PxS P 2 

n 7 


+ P | 2 ... ^. .. b . + P 22 J^L- 

Sx 2 d P | Sx 2 Sz 

+ P 2 ^ ■^■ 3 K ... + 2 P 2 + P | + 2 Pl 2 /" 

3z3xi 3x 2 Sz 2 3 z 3 \5z5x 2 < 

+ 2 P 22 (■ + P 2 ■' — + P 12 

\SzSx 2 5pa Sz 2 ^pa J SzSp^ 

+ P | 2 J3L + P22 &u 

a z a P i az 2 

r J / 

^ + 2 P s ^ + pi ^ + 2p l2 ( . ^ 3 W + p 2 ^ 

a Pl ax| a Pl ax 2 a z a Pl az 2 \ax 2 ap 2 a z a P 2 y 

2 P22 f g .. - + v s - ^ a + pfs + “ 

\dp 1 OX 20 p 2 opx^zr 


P 12 


> ( + p 2 c )3 W 

\a z ax 2 a Pl a z 2 a Ply 

-^S_ + 2p 12 p 22 _3fw 

a z a P f a z a Pl a P 2 


i. 


+ PI 2 


+ 2 Pi2P 22 ^ 
a P f a P fa P 2 


zS Pl 


a^ x „ a% 

^pi^pI s Pl a z 

+ P22 f-^L_ + 2 + 2 P2 + 2 p 2 + p| J!w_ 

dp2^ x 2 dx^dz Sp^Sx^Sz dz 2 Spa^z 2 

+ vy 1= f 83w * ^ *Vs ^ 

\a P 2ax 2 a Pl a z a Pl ap 2 a z a; 

♦ 2 P 22 f . a °» . ■■ + J&l + P 2 -i*-) + pf 2 

\ax 2 a P i a z a P 2 a z a P 1 / a P 2 a] 

+ 2 Pl2 P22- + P§2 + P22 ^ - 

^Pi^P2 a P f a P 2 a s 

r 

+ 2 Pi2 2 'j — ^ + P2 — + Pl2 

dx^dpx ozopx ^Pi 


d W 

a P2 ap! 


^ w 
Spa^z 

+ P 22 ..- 5 -W- 

a P f a Pl a P 2 
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+ P 222 


+ P 22] 


+ P222 


s( 

- + P2 

a 2 w 

- ] + 2 Pl2 

a 2 w 

Vdx 2 a P£ 

a z a Pi 

V 

a Pl a P2 

a 2 w A 

a 2 w 

"Ho 4* 

S 2 w + 

a 2 w 

ax 2 a Pl 

^zS Pl 


dpf Sl2 

ap 2 a P i 

r- 

a 2 w + 

a 2 w 

P2 + 

a 2 " Pis 

, + if" 

ax 2 a P2 

aza P2 

a Pl a P2 

a P f 


, o„ 5 2 w , aw 

+ 2 p 22 — 5 - + 

OP 2 5z 


P22 


(4.220) 


Wow substituting Eqs. (4,217) -(4.220) into Eq.. (4.216) and eliminating 
P 222 'by Eq. (4.210), we get 


f O + f lPl2 + f2Pl2 + f 3 P22 + f“4PllPl2 + f 5 P§2 
+ fePl2.P22 + frPl2P22 + fsPl2Pl2 + fsP22 

+ floPl22 + fllPll + f l2PllP22 + f 13Pl2 

+ f" 14 P 22 1 = 0 


(4.221) 


where 


:■ <*L + P2 + j| a? 2 _ . p, . p,p 2 

dpi c3xc)y dxdty Sy 2 3y3b 


a a w ~ 2 S 3 W 


\3 X T 7\3 T 

Sy 3 3P2 


3Pi 


dydf 1 


P2 f pf ^ + P2 

V a^y 2 


, 0^2 a 3 W , . j (, d s W , B 2 W , 

+ + ^(5 — + ^ — - 


SySps 


ata P2 at. 


(4.222) 


- 2 p 2 


a a w 

apia P2 


_ a 2 w , 

Pi y-g “ 3 

OPl 


a 3 w , ■ a 3 w 

+ P2 


a y a P f 


at apt 


(4.223) 
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P2 aw + aw _ a P2 fJ fw. + P2 -ifW- 

St Sy l Sp 2 Sy dp 2 St 


_ 4, S 2 ¥ . ^ S 2 ¥ 

1 \ Sy^Pi 2 StS p ly 


2 ^^— + -n ^ 

ySy^Sp! 2 SyttySpi 


S 3 ¥ 

Bpi^y 2 


SP2 f .a 3 * + P2 ^¥ . . 

VStSySpx Sf 2 Sp! 


^PiSySt 


pi a3w + 5^f a 2 w 

dpiSt 2 SpiBps 


(4.224) 


S g ¥ 

SxSp 2 


4 ? WVY , __ W VV , „ / W vv , 

f 3 = n— + P2 + Pi + 2 P2 


S 2 ¥ 

SySp 2 


S g ¥ 

StSp 2 


P ( S 3 ¥ .,. „ S 3 ¥ 

V Sy 2 Sp 2 2 SyStSp 2 


S 2 ¥ f S 3 ¥ S 3 ¥ \ cfw S 3 ¥ 2 S 2 ¥ 

SySt ^ 2 \ StSySp 2 P2 St 2 Sp 2 / P2 St 2 Sp^y 2 SySt 


+ gp2 a 3 w + gpe gg t pi a 3 w _ 

Sp2^ySt St); 2 Sp 2 S\|; 2 


+ Pl ™ + 3^ 9* 
St Spf 


(4.225) 


(4.226) 


f = 2 — - 5 f S 3 ¥ + S 2 ¥ + S 3 ¥ \ 

Pl Spf l SySpl StSp 2 P2 SfSp! y 


(4.227) 


-2p 2 £* + Pl -^L_ 

Spf s Pl s P£ 


/ _^w_ + P2 a 2 ¥ 

v Sp 1 SySp 2 S Pl StSp 2 


u w 

SpiSt (4.228) 


SpfSp 2 


(4.229) 


SpiSpl 


(4.230) 


(4.231) 
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£ 10 - 1 


*■» + 2p2 

8 2 ¥ 

6y8p x 


^ 11 = P2 ^ 

^ S 2 ¥ 4 . 


s.dpi<6y 

f 12 

= P2 : 


1- gp sa y w > 

5pf Bpi^paJ 


P2 


dPity, 


- 13 “ 


6 2 ¥ 

^PlSp 2 

6 3 ¥ 

6 P f 


(4.232) 


(4.233) 


(4.234) 


(4.235) 


-14 


6 2 ¥ 

BySpx 


8 2 ¥ , d 2 ¥ _ , 8 2 ¥ „ 

P2 + — Pl2 + P22 

^j/dp! opf 5 pi^P2 


(4.236) 


Since ¥ is a function of x, y, \|r, p x , and p 2 , Eq. (4.221) is satisfied 
identically only if all the coefficients are zero. Therefore f 0 ,...,f 14 are 
all equal to zero, which gives fourteen equations 


f O — f 1 — ooo — f 14. — 0 


14 


(4.237) 


from which the functional form of ¥ is determined. Prom Eqs. (4.234) ¥ can 
be separated into two terms as 


W(x,y,\|f>p 1 ,p 2 ) = W 1 (x,y,f,p 1 ) + ¥ 2 (x,y,\|r,p 2 ) . (4.238) 

Equation (4.226) implies ¥ x is linearly dependent on p x . Therefore, 

Wi(x,y,f,p x ) = W xx (x,y,\|r)p x + W 12 (x,y,t) . (4.239) 

Equations (4.232), (4.233) and (4.236) then give the same relation as 
follows : 

+ p 6¥^ = 0 (4.240) 
oy * 


94 



Since W 13 _ is independent of p 2 , it means ¥ 13 _ is independent of both y and ijr 
and this is a function of x alone. Equation (4.223) is seen to be satisfied 
identically. ¥ now becomes 


¥ = ¥ 11 (x)p 1 + ¥ l2 (x,y,t) + W 2 (x,y,i|r,p 2 ). (4.24l) 

Substituting this form of ¥ into Eq. (4,227 ) , we get 

2 Pl 3fik . 3 f. ' + + p 2 J^) , 0 • (4.242) 

dp! \ dy 2 dpi d\|/dp 2 dtydpf J 

Since ¥ 2 is independent of p x , Eq, (4,242) gives 

^ = 0 (4,243) 

op 2 

and thus 

^2- = 0 , (4.244) 

Equation (4,243) implies W 2 is linearly dependent on p 2 o Also, from Eq, 
(4 0 244), the coefficient of p 2 is independent of \|r. The characteristic func- 
tion^ W, can he written as 

W(x,y,\|/,Pi,P 2 ) = W 11 (x)p 1 + ¥ 2l (x,y)p 2 + ¥ 3 (x,y,t) . (4.245) 

Substituting this form of ¥ into Eq. (4.244) we get 

^3. + p 2 _ dW 2j . + d ¥ 11 \ _ o e (4.246) 

dy \d\|r dy dx J 

Since ¥ 3 , and ¥ 2l are all independent of p 2 , Eq. (4,246) is separated into 
two equations, namely, 

^3 = 0 (4.247) 

dy 

_ kki + dW n = o . (4.248) 

dijr dy dx 
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Equation (4.247) shows that ¥ 3 is independent of y. Since hoth ¥21 and ¥ X1 
are independent of f, Eq. (4.248) means ¥3 is linearly dependent on i|r. There- 
fore, 


¥ = ¥ 11 (x)p 1 + ¥ 2 x(x,y)p 2 + ¥ 31 (x)f + ¥ 32 (x) . (4.249) 

Employing this result in Eq. (4.248 ), we get 

+ ¥ 31 = ffi& IL . (4.250) 

dx 8y 

The left-hand side of Eq. (4.250) is a function of x alone , which in turn 
means ¥21 is linearly dependent on y, i.e., 


w 2 i(x,y) = W 2ll (x)y + ¥ 2 12 ( x ) . (4.251) 


Equation (4.250) then hecomes 


d¥ 


dx 


H + ¥ 3 ^ — ¥2 


11 


0 


(4.252) 


which will he used later. 

The characteristic function now hecomes 


¥(x,y,i|r ,Pi,p 2 ) = ¥ 11 (x)p 1 + [W 2ll (x)y+W 2l2 (x) ]p 2 + ¥ 31 (x)\|r + ¥ 32 (x) . 

(4.253) 


This form of ¥ will satisfy Eqs. (4.228)-(4.231 ) and (4.235) identically. 

¥e have two equations left, namely, Eqs. (4.222) and (4.225). Substituting 
¥ from Eq. (4.253)' into Eq. (4.222) gives 

fi'Wn + pl ^ + p| d ~ W31 + ( J¥ 2 n+¥ 3 2. ) = 0 (4.254) 

dx dx 

which can, in turn, he separated into two equations: 
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dx 


0 


(4.255) 


+ 


d¥3i _ 
dx 


+ ^(3W 211 +w 31 ) = o . 

Similarly, putting W from Eq. (4.255) into Eq. (4.225) gives 


(4.256) 


K pi + 

<i W 211 y + dw 2l2 

P2 + dW 31 

\|r + dW 32 

l dx 

_ dx dx _ 

dx 

dx 


+ p 2 | d ^ 2 ii y 
\ dx 


dWg 2 . 2 . 
dx 


+ Pi(W 2ll -W 3 i) - 0 


(4.257) 


which then is separated into 


dW 1T 

dx 


dW 31 = 0 

(4.258) 

dx 

dW 3 2 __ n 

(4.259) 

dx 

W211 + W31 “ 0 0 

(4.260) 


Equations (4.258) and (4.259) show that W 31 and W 32 are constants and so from 
Eq. (4.2.55), W 211 is also a constant. Equations (4.252) and (4.260) give the 
same result, namely, W u being linearly dependent on x. 

To conclude, the characteristic function W can be written as 


W(x,y,i|f,Pi,B2) 

= (ao+axx)?! + [g(x)+b3_y]p 2 - [co+cxtl 


(4.261) 


where a Q , a 1# bj, c G , and c x are constants and g(x) is an arbitrary function 
of x. One relation relating <J> [defined in Eq. (4.221)] and the constants 
can be obtained from Eq. (4.256) as 


(ao+aix)^ 1 - (c 1 ~5b 1 )5z( = 0 


(4.262) 
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and also, Eg.. ( 4.260 ) furnishes another relation among the constants, which is 


a i - t>i - c x = 0 . 


(4.263) 


The absolute invariants and the restriction on f will he discussed in 
the next section. 


4.3.2 Limitation on the Mainstream Velocity 

In the introduction the point was made that the type of problem involves 
the determination of an arbitrary function, <f>. For such cases, <not only the 
characteristic function ¥ must be of the form given by Eq, (4.26l) and Eq, 
( 4 . 263 ) satisfied, but the arbitrary function 4 has ho satisfy Eq., (4.262) 
if the given differential equation is to be invariant. 

From Eq.. (4.262 ) , 


(ao+a^x)^’ - {c x -yo x )<f> = 0 . 

Upon integration, for a x ^ 0, a 0 =0 

4 = (constant )x^ Cl_ ^"' :)l ^ / ^ C;L+ ^ ;L ^ (4.264) 

From Eq. (4.211), 

1^1 = (constant )x ( c i"3hi )/( c i + ^>i ) ( 4 . 265 ) 

2 dx 

we then get 

u e = (constant )x^ Cl ”"* Dl ^^ Cl+ ^ 1 ^ (4.266) 


This is seen to be the mainstream velocity associated earlier with the linear 
group . 


For a x =0, a D / 0, 


c i~3bi 

(constant )e 0 


(4.267) 
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and thus, using Eg. (4.211), 


c 1-3*1. . 

u e (x) = ( constant )e ^ a o 


(4.268) 


which is seen to he the mainstream for the spiral group. 

For the general case in which both a Q and a x are nonzero, the mainstream 
velocity is found to be 


u e (x) = (constant ) (a^+ao )^ Cl+l31 )Aci+b x ) (4.269) 

V- 


which again belongs to the linear group. 

We therefore conclude that the mainstream velocity, u e (x), belongs either 
to the linear or the spiral group, i.e., powers or exponentials of x. Wo 
other forms are possible —a point long assumed but not proved. 


4.3.3 Absolute Invariants and the Transformed Differential Equation 

With the characteristic function, W, obtained as in Eg. (4.26l) we now 
make use of the general theory to find the absolute invariants and the trans- 
formed differential equation. We shall expect the invariants to correspond to 
linear and spiral groups. From Eg. (4.265), we obtain 

^ (4.270) 

Oh a 2 m 

where the functions a x , Cfe, and m can be obtained by putting Eg. (4.261 ) into 
Eg. (2.13) which then gives 

= - T -^ = . (4.271) 

a 0 +a x x g(x)+b!y c 0 +c x ty 

Since this equation has two independent solution^, we have two absolute 
invariants. Let us consider two special cases. 


Case 1. a Q = 0, g(x) =0, c 0 = 0 
Equation (4.271) becomes 
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dx 

a x x 


(4.272) 


dy _ d\j/_ 

t>i y c x t 

The two independent solutions to Eq.. (4.272) are 


and 


x t >i/ a i 


constant 


c i/ a i 

x 


constant . 


The houndary-layer equation can he expressed in terms of these two invariants, 
Thus^ we get the similarity transformation 


n = — and f(n) = — k — • 

x ^x/ a i x c i/ a i 


(4.273) 


The mainstream velocity for this case is given in Eq.. (4.266) as 


U 6 (X) = kjFi-h^o+h;,) 


(4.274) 


where k x is a constant. Recalling that Eq.. (4.263) has to he satisfied, Eqs. 
(4.273) and (4.274) can he written as 


” = and f(n) = 

and 


u e (x) = k^ 111 


(4.275) 


(4.276) 


which is seen to he in the same form as the so-called Falkner and Skan sim- 
ilarity transformation (see Ref* 6) 0 The transformed differential equation 
is well-known and can he written as 


f"’ + | (mfl)ff" + m(l-f’ 2 ) = 0 . 


(4.277) 
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Case 2. a x = 0, g(x) = 0, c Q = 0 
Equation (4.271) becomes 


dx _ fly, _ djjr 

a o _ iiy ~ c x t 


(4,278) 


The two absolute invariants can he found to he 


- ia 


and f ( rj ) = 


t 


(4.279) 


The mainstream velocity is given hy Eq., (4,268) as 


u e (x) = k 2 e^ X , 


(4,280) 


The transformed differential equation is 


f n» = p( f «a _ 1 ff ") _ pk | . ( 4 . 281 ) 

2 


The transformation is seen to he the spiral group. 

The more general cases in which hoth a 0 and c 0 are nonzero in Case 1 
and c 0 is not zero in Case 2 poses no problem. These assumptions merely in- 
troduce an extra constant in x and in the transformation [cf, Eqs. (4,266) 
and (4,269)]. However, the case in which g(x) is not zero needs further in- 
vestigation. For this case, Eq., (4.271) becomes: 


dx 

a 0 +a!X 


dy _ d\[r 
h x y+g(x) Cit+Co 


(4.282) 


Following the same steps as in Case 1, the absolute invariants are found 

to be 




y 

(a 0 +axx) ai 



g(x)dx 

+ 1 

(ao+axx ) 8,1 


(4.285) 
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and 


f(Tl) 


c it+Co . 

(ao+a ^) 01 &1 


The mainstream velocity is given by 


u e (x) 


,(a 0+ a lX f 0l - b ^° 1+ ^ 


based on Eq. (4.2 69 ). 

The boundary-layer equation is transformed to an ordinary 
equation as 

Cl " ' b - 3 * f ,2 - ff" = f" ! + k§ » 

a i a x c i+*i 

By setting 

u e (x) = ksCao+aix) 111 


then 

m = . 

ci+bi 

Also, Eq„ (4.2 63 ) has to be satisfied, i.e„, 

a i = "bi + Cx . 


Equation ( 4 . 283 ) bo (4.286) become 


y 

r g(x)dx 

(a 0 +a!xj 

F#) J 

(a Q +axxj 

p-m 


f(Tl) 


C 1** C 0 _ 

, .l+m/2 

(a 0 +a x x) 


(4.284) 


(4.285) 


differential 


(4.286) 


(4.287) 


(4.288) 


(4.289) 


(4.290) 


(4.291) 
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and 


f"* + mkf - mf ' 2 + ^ ff" = 0 . 


(4.292) 


For g(x) = 0, this is seen to reduce to the Falkner-Skan's flow. In order to 
transform the boundary-layer equation to an ordinary differential equation, 
g(x) can be any function of x; however, for g(x) £ 0, the boundary condition 
at y = 0 cannot be transformed. We therefore conclude that for incompressible, 
two-dimensional, laminar boundary-layer equations, the linear and spiral groups 
are the only two groups possible. The group represented by Eq. (4.282) will 
transform the boundary layer equation from partial differential equation to 
an ordinary differential equation, it fails, however, to satisfy the boundary 
condition at y = 0. 


4.4 SIMILARITY ANALYSIS OF THE HELMHOLTZ EQUATION IN GENERAL CURVILINEAR 

COORDINATES 

In the present article the role of a coordinate frame relation on a sim- 
ilarity analysis is considered by examining the two-dimensional Helmholtz equa- 
tion in general curvilinear coordinates. While the present method of analysis 
applies to both linear and nonlinear partial differential equations, there are 
at least two reasons for choosing this particular equation for presentation 
of the present method of analysis. The first is its simplicity and its wide 
use in the potential theories of a large number of physical problems. The 
second, and the most important, reason is that the separability of this equa- 
tion has been studied thoroughly by Moon and Spencer. 43 By separation of 
variables we mean the original partial differential equation is reduced to 
two ordinary differential equations each having one of the original indepen- 
dent variables as an independent variable. In the case of similarity, the 
original partial differential equation is reduced to one ordinary differential 
equation where the independent variable is a function of the original inde- 
pendent variables. The present analysis is a parallel study to that of Moon 
and Spencer43 on the conditions of separability. Here we want to derive the 
conditions under which similarity is possible by requiring that the Helmholtz 
equation be invariant under the infinitesimal transformation. The basic prob- 
lem is the determination of the characteristic function, W. The difference 
between this example and earlier examples is that the resulting differential 
equations for the solution of W involve unknown metric components of the 
curvilinear coordinates. These equations form the conditions for the exist- 
ence of similarity solutions. They can be used in two ways. In the first 
case where a coordinate system is given, similarity solution is said to exist 
if substitution of the known metric components in these conditions result in 
solutions for the characteristic function, W. Once W is known, the searching 
of all possible groups of transformation can be made by following the same 
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steps as in previous examples. In the second case where a group of transforma- 
tion is given for which the characteristic function, W, is known, these condi- 
tions give limitation on the metric components for the general curvilinear co- 
ordinates.. If substitution of the metric components of a given coordinate 
system satisfies these conditions, then similarity solutions exist for the co- 
ordinate system for that group of transformation. In both cases, of course, 
the boundary conditions have to be transformable for true similarity to exist. 


4.4,1 Review of Separability Conditions Discussed by Moon and Spencer 

The following is a short review of the separability conditions discussed 
by Moon and Spencer. 13 This is done for comparasion with results to be ob- 
tained from the theory of continuous transformations. We consider two dimen- 
sional form of the Helmholtz equation given by 


id/ g 1 / 2 dg$ \ | d / gi/ g 5^ \ 
g 1 / 2 dx x \gxi d dx 2 V^g 2 a dx£> J 


+ OL^i 


0 


(4.293) 


where g n and g l2 are components of the metric tensor and g is the determinant 
with elements g^ . 

Now, let the unknown function be expressed as a product of two functions: 


<f> = UiUiJUsCxe) . 


( 4 . 294 ) 


Substitution of Eq. (4.294) in Eq. (4.293) gives 

. 1 d 


V 2 


1 d f g 1 / 2 dU 1 ' 

u i dx x Vgn dx ly 


U 2 dxsVgsa 



+ a x = 0 . (4.295) 


The most general condition that will allow separability is that g 1 / 2 /gj_p 
is a produce of t*wo functions; 

si/i = f 1 (x i )P 1 (x 2 ) (4.296) 

gll 
g 1 / 2 


f 2 (x 2 )p 2 (x 1 ) 


( 4 . 297 ) 



Substitution of Eqs . (it-, 296) and (it-. 297) into Eq. (it-. 295) gives 


1 d 

gi/ 2 U i dx i 



+ SB. -A. L 

U 2 dx 2 y dx 2 y 


+ a x 


0 . 


(4.298) 


In this equation, g, F± and f j_ are determined by the coordinate system and 
are entirely independent of the boundary condition that characterize a partic- 
ular problem. But in the solution, (f> = U X U 2 , the U's are functions of both 
the coordinate system and the separation constants. Now, suppose we dif- 
ferentiate Eq. (4.298) with respect to Qd and a 2 , we get 


1 Sod 


i_ _1_ f f x 

Ui dx x V dx ly 


f+ Fs: 


all d 


f2 


dU s 


8ax U 2 dx 2 \ dx 2 


+ *g 


1 / 2 = 


(4.299) 




1 d 


dU n 


Ui <3- x i \ X cLx ly 


+ Fg 


8o!2 


1. Es. 

U 2 dx 2 V dx 2 


= 0 


(4.500) 


Introducing the notation 


^ij( x i) = 


1 a 


fi( Xi ) a XJ 


_i (f . Ei 

Ui dXi V 1 dxi 


(4.501) 


Equations (4.299) and (4.500) becomes 


flFl^nUx) + f2F 2 ^2l( x 2) = g 1 / 2 ( 4.502) 

fxFxSz(x2( x i) + f‘2F 2 jz( 22 (x 2 ) = 0 . (4.505) 


Equations (4.502) and (4.505) are solved for fxFx and f 2 F 2 and the results 
are (if s / 0) 


fxFx = g 1 / 2 ^22. 

s 


(4,504) 


f 2 F 2 = - g 1 / 2 iis. 

s 


(4.505) 
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•where 


^11 

5^12 

4 2 i 

9^22 


(4.306) 


which is known as the Stackel determinant. 

Comparison of Eqs. (4.304) and (4.305) with Eqs. (4.296) and (4.297) shows 

that 


gn 


s 

5^22 


(4.307) 



(4.308) 


Equations (4.307) and (4.308) are called the first condition for simply 
separability. 


Also^ from Eqs. (4.304) and (4.305), 


g x / a 

s 


flUx) 


ps * 

F l (x 2 ) 

= f o ( Xo ) 

r , , n 

Fi(xi) 

4 22 ( x 2 ) 

•‘-id \ ed J * 

-5z! 12 (x 1 ) 


This is possible only if 


Si/f. = f 1 (x 1 )”f 2 (x 2 ) . 


(4.309) 


(4.310) 


Equation (4.310) is called the second condition for simple separability. 

It can be shown that the two conditions are both necessary and sufficient 
for the simple separability of Helmholz’s equation. Detail of the proof is 
given in the original work of Moon and Spencer. ^3 

We now consider an example in which the cylindrical coordinate is con- 
sidered, for which 


x = r cos Q, y = r sin 0 


(4.311) 
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From Eqa, (4,307) and (4,308), 


= 1 


and 


s 

9^12 


= -r 


If ^22 i s taken to toe 1, then 

s = 1 and <j > X 2 

However, s is defined in Ecu (4.306) as 


s = 


^11 5^12, 


9^ii 

1 

“ r 2 

S^2l 9^22 


9^21 

1 


For s = 1, one possibility is 

<j> xx = 1 and 52(21 - 0 


The second condition then heoomes 

st/1 

s 

This condition is satisfied if 


= r = fi(r)f 2 (0) 


f i(r ) = r and : f 2 ( 6 ) = 1 


We therefore conclude that the Hel m holtz equation in 
ordinates is separable. 


(4.312) 

(4.313) 

V. 

(4.3l4a,to) 

(4.315) 

(4.316) 

(4.317) 
(4. 318 a, to) 

cylindrical co- 
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4,4.2 Similarity Analysis of the Equation 


Let us put Eq. (4,293) in a slightly different form as 


Jl( ^ + Z(h i 

OX! V h^ 0X3/ 0X2 \ t L2 0X2 


It - ) + aitete/ = 0 

0X2 / 


(4„319) 


where the following relations have teen used: 


l/2 = 


kite, S 11 = kf. 


= hf 


¥e now want to find the conditions on the two unknown functions h x and 
h 2 which make similarity possible. In other words, the conditions thus ob- 
tained will enable one to decide If similarity solutions exist for the co- 
ordinate system under consideration. 

Carrying out the differentiation and rearranging the terms, Eq. (4,319) 
becomes 


fy, v, 2 ^h 2 , o 

ihlh2 . ^ 


,|4 Pl+ ( hShs g, 


v3 dh 2 , ^ 
“1 ) P2 




+ hihfpn + hfh 2 p 22 + Q!h,fh|^ = 0 


(4.320) 


where 


Pi 


88 88 8 2 8 
= £? P2 ■ sb Pl1 = aJ’ P22 - * 


8 2 ^ 

8xf 


(4.321) 


We now make the infinitesimal transformation: 


*{ = Xi + e| i (x 1 ,x 2 ,^,p 1 ,p 2 ) 

4 X - 4 + e^(x 1 ,x 2 ,^,pi,p 2 ) 

P 1 = p i + e5t i( x i>te,^,Pi,P2) 

Pjj = P ± j + €rtij( x i,Xs,^Pi,p 2 ,p 1;L ,Pi 2 ,p 22 ) 

(i,J = 1,2) (4.322) 
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where the transformation functions and rc 1 j can he expressed in terms 

of a characteristic function, W, as 


e - &W 

El “ 55T 

C - p, - W 

OPi 


-«i = | 2 - + p, * 


1 ^ 


(4.323a) 

(4.323t>) 

(4.323c) 


- nil = * 2 Pl -J5!_ + pf + 2 Pll (2^W_ + Pl J50\ 

a x f 1 a Xl a$z( 1 a ^ 2 11 la Pl a Xl a^a Pl y 


+ 2 Pl2 h-Sfw. + Pl Jfw_\ + jf, ^ + 2 PllPl2 Jfw_ 

ia Xl ap 2 a P2 <v y a P f a Pl a P2 


+ P f 2 a?w 

a P f a^ 


(4.323d) 


a 2 w 


a 2 w 


-rt 22 = — — — + 2p 2 — — — — + p 2 ^ ^ + 2 Pl2 ( — — — - — + p 2 

x^ x^2 1 1 ' " * 


A 


a 2 w 


a 2 w 


a X2 a^ a^ 


V Sx 2 Spi d^dp 3 


+ 2 P 2 2 1 r— ^ 7 — + P 2 + Pl 2 t 2p 12 p 22 


a X 2 a P2 a^a P; 


Spf 


aw 

s Pl a P2 


2 a 2 w . _ aw 

+ ppp — — + ppp ■ — ■ 

dpi ^ 


(4.323e) 


The Helmholtz equation, Eq. (4.320), denoted hy F = 0, is invariant under 
the infinitesimal transformation if 


UF = 0 . 


(4.324) 


Expanding the operator in full^ we get 




aF 


+ 


aF 


1 a P . ^ a P 


= 0 


ij 


(4.325) 
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Patting the I^elmholtz equation, Eq, (4.320) into this expression, we 
then get 


(!r) 

2h x h 2 + h f Shs dhi + hfh 2 

Sx x Sx 2 Sx x Sx 2 dx x ^X 2 


1 Shg. <^1. - h| S 2 hi 
Sx x Sx x Sx§ 


Pi 


,2 Shi Sh2 _ h? _^2 


- 3hl Piii 


3 Sh x ,, ,2 Sh 2 

12 st + 3hilte sr 


Pll + 


3h!h| gl ♦ 3hHJ || 


Shi + h 3 

OXJ. dx, 


Sx x Sx 2 Sx x Sx 2 


P22 


P2 


+ 12 { 




^2 Sh n Sh 2 , 01 ^ v Sh 2 Sh 2 , .z S 2 h 2 

h 2 — - ^ — — + 2 h x h 2 r— t -r — + hih 2 r — r — 

OX^ OXjl OXjp OX^OX^ 


2 Sh^ Sh^ _ h3 S g h x 


3h 2 
2h x h 2 


Pi 


bx.2 Sx x “ S2feSx x 

dhi\ 2 _,_ r,2 Sh 2 Shi , v2r, s 2 h x 2 Shj^ Sh 2 
EiiJ, ) + hi sr-s- + h x h 2 - 3h x ^ - — 


SXj3/ 




Sx 2 Sx 2 ' i+Xi ^ Sxl Sx 2 Sx 2 

/ 2 , Sh x 3 Sh 2 \ 

il2 + hi s^) p “ 

+ £ah^hf 


- h' 


3 S 2 h 2 
1 Sx§ 


P2 


») Pl1 + ( ! 

+ <* ( 3h M |j£ + 3^1 |j) 

+ ** ( hihS i: " hl it) + ” 2 ( h!te i& " hp 

+ it xx h x h 2 + rt 22 h x h 2 = 0 „ 

Putting Eqs, (4,323) into Eq.. (4.322) and we get: 

f 0 Pll + flPl2 + f2Pfl + f 3P 1 lP 12 + f 4Pl2 

+ f 5 P 12 P 22 + f 6 P 22 + f > 7 P 22 + fa = 0 


(4.326) 


(4.327) 


where 
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dh n , n _2 5ll2 \ , SW 


^ v 112 st * 3hihS srr 3 p. 


h x h| 


n— l + 3hih 2 ^ — 
dxs ox 2 


2 ( ^ 

a 2 ¥ N 

1 + ^ 

_ \apia Xl 

Pl a^a P J 

1 -+• — 
c)^ 


■1 


hi hi 


a 2 w 


+ pi 


a 2 ¥ 


, c)Xidp2 ^P2^j^ 


h?h 2 


2 /1£w_ +P2 JltL 

VBxidpx dyfopi 


f - h h 3 S 2 W 

12 — “ -h-l-h-2 N 2. 

opi 


-3 - 


-2h a h| 


l 4 - 


-hihg 


a 2 w 

a Pl a P2 

hfh 2 


3 a 2 w v,3v, a 2 ¥ 


Spi Spf 


-2hfh 2 ^ 


dpi&Pa 


L 6 


-hfhs 


a 2 ¥ 

apl 


-7 - 


aw_ 

Spi 


3hfh 2 


ah! 

a Xn 


.3 ah 2 \ , a¥ f - 7 i^ 2 v,_ ah x + ^3 ah 2 


+ — ( 3h 2 h 2 ^ hi 
SP2 \ 


+ h X ~g. . T I J11J.U3 - T u x 

ax n j a-D? v 0x2 ax 2 


3 f a 2 ¥ a 2 ¥ 

2h!ha + P£ 


;) - | 


(4.528) 

I 

(4.329) 

(4.330) 
(4.35i> 

(4.532) 

(4.333) 

(4.334) 

(4.335) 


ill 



U^jL^x OV, v_ a. V ^2 f¥he 




*1,2 Sh 2 Sh x t. 3 S 2 h T _ 

5112 3*1 st ' 112 SSf Pl 


5hi SliT , ,2 c)h s ^h.T , v2r a 2 h-| 
Sx 2 1 3x x Sx2 1 2 SXi^Xs 


, h s dh T dhg _ h 3 a hg 
1 Sx x 5xg 1 dx-ictecp 2 


+ ctf 3hf hi + 3 hf h| 


3-u2 ^hg 

lh2 SSI 


h g ii S-a. + 2^ ^ 3i2. + h lh | ;iia. 

OXg OXj 


Sxg Sx 3 


9^2 3x x 


,., 2 Shp ShT -, 3 d hi 
3hg 5-1- - hg - — p x 

oxg oxi ax^oxx 


+ 211^2 




^2 Shn 3hg v 3 S 2 hg ,, 
5hi — “ r-®- “ h x ~-^- pg 
oxg oxg a;xg 


+ cuf ( 3hfhl + 3hfh| & 


+ Pl §~ + P2 §~- - W Qhfh| 


5 2 W . g S 2 W 


hl hi / — + 2 Pl - 


,3, f S 2 W 0 S 2 W . 2 S 2 W 

12 Uxf 2ps Sx £ ^ P2 ^ 2 . 


aw . aw\ /T .2 ahg , 3 ^ 

Sxi + ^s?) ( h ^ ssr - ^ sj; 


— + P2 — hf hg ^ 

axg aa/ V ^ x s 


1 a X g 


in Eq_ 0 (4„327) can he eliminated, hy using Eq.„ (4.320) which, !] 


shorthand form, can he written as 


‘(gn+SlPll) 


(^.537) 


P22 
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where 


So 






p 2 + Ohfhl^ 


K^.338) 



(4.339) 


By substituting Eq. (4.537) into Eq. ( 4 . 327 ), we get 


HqPh + HiPis + Ifepfi + H3P12 + H 4 = 0 ( 4.340) 

where 


H 0 

= f Q + 2 f s g 0 g 1 - fygi 

(4.341) 


Hi = fi - f 5 g 0 

(4.342) 


H 2 = f2 + feSi 

(4.343) 


h 3 = f 4 

(4.344) 

H4 

= f 8 Sq^Y + ^ 6 g§ » 

(4.345) 

Since the characteristic 

function, W, is a function of x 3 , x 2 , 

<f>, Px, and 


p 2 , Eq. (4.340) is true only if the coefficients of each term are zero re- 
spectively; that is, 


H q = h x = H 2 = H 3 = H 4 = 0 o 


(4.346) 


The set of equations, represented by Eq. (4.346), forms the conditions 
for the existence of similarity solutions. These conditions can be used in 
two ways. In the first case where a coordinate system is given and we want to 
know if similarity solutions exist or not. Since now h x and h 2 are known, 
they can be substituted into these conditions and if there exists a solution 
for the characteristic function, ¥, similarity solutions then exist. The 
searching for all possible groups can be made by following exactly the same 
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technique as given in the previous examples. In other words, as long, as sub- 
stitution of h x and h 2 for a given coordinate system into these conditions 
results in solutions to the characteristic function, ¥, the coordinate system 
under consideration possess similarity solutions. 

In the second case, we consider a given group of transformation and we 
want to know the conditions h x and h .2 must satisfy for similarity to exist. 

For a given group of transformation, the characteristic function, W, is known. 
Substitution of this form of W into the conditions, Eq. (4.346) will result 
in five equations connecting hi and h 2 . They are the conditions that a given 
coordinate system must satisfy if similarity solutions exist for this given 
group of transformation. 

Two examples are given below. v 


4.4.3 Conditions of Similarity for a Given Coordinate System 

As a simple example, consider the rectangular coordinate system in which 
both hi and I 12 are unity. Then the conditions for the existence of similarity 
solutions, Eq. (4.346), give 




+ pi 


5 w 

8^8p n 


8 ¥ 


a<f + 


f 8 2 ¥ 


dpi r V dx 2 dp 2 


+ P2 


8 ¥ 

cSj^Ps 


= 0 


(MV 7 ) 


8 g ¥ 


+ Pi 


8 g ¥ 

dp 2 d^ 


8 a ¥ 

. dx 2 dp 3 


+ P 2 


8 a ¥ 

6^p 3 


+ of 


8 g ¥ 

dpidp 2 


= 0 


8 2 ¥ , d g ¥ 

8p? dpi 


= 0 


(4.348) 

(4.349) 


8 ¥ 


8 ¥ 


a|Pl diT +Pg 

' 8 2 ¥ 


x t » , d 2 ¥ , - 8 2 ¥ P 8 2 ¥' 

Wl ST 2pi Pl 3? 


+ 2p 2 -ifSL + pi 

8xf 8x285^ 8$z( 


8x2 dp2 d^dp2 


8j? 


8 2 ¥ 


oF( 2 *-£• = 0 . 
dp2 


(4.350) 
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This means if solution to Eq_s. (4. 347) -(4. 350) for W exists, similarity solu- 
tions will then exist for this given coordinate frame. The searching for all 
possible groups of transformation satisfying these conditions shall not he our 
concern. The only conclusion being sought is whether similarity solutions 
exist and if they do what are the conditions the characteristic function, W, 
must satisfy. Equations (4.347)- (4.350) form those conditions and by trying 
various functions for solutions, the conclusion is easily drawn that simi- 
larity transformation do exist. An example will be given below. 

Consider now the following form of the characteristic function, ¥, 


W(x 1 ,X 2 , J z(,p 1 ,p 2 ) 

= W 1 (x 1 )p 1 + W 2 (x 2 )p 2 + w 3 (x 1 ,x 2 ,^) . (4.351) 


For this form of W, Eqs. (4.348) and (4.349) are satisfied identically and 
Eqs. (4.347) and (4.350) give 



Equation (4.352) gives 

Wi = Ci*i + c 2 
W 2 = c x x £ + c 4 . 

These are substituted into Eq. (4.353) and we get 

J 0 Pi + JiP£ + J 2 P 1 + J3Pi + J 4 = 0 


(4.352) 


(4.353) 


(4.354) 

(4.355) 


(4.356) 
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where 


J x = 


J 2 = 


J 3 = 


^3 

3_ 

3 x 2 3 jz ( 

SfWa 

8 2 ¥ 3 

d^ s 


J 4 — -C£W 2 


3 2 ¥. 


3xf 


a - 


^Ha + 2ctf Cl + cct^f 
3x1 op 


(4 '.357a). 
(4.35Tb) 
(4.357c) 
(4.357d) 
(4.357e) 


Since all the J's are independent of the p’s, Eq. (4.356) is satisfied if all 
the J’s are zero; that is, 


Jq — Ji - J 2 - J 3 - J 4 _ 0 • 


(4.358) 


The first four conditions in Eq. (4.358) indicates that the function 
¥ 3 , should he of the form 


W3 - + W 3 i(x 1 ,X2), 


(4.359) 


Substitution of Eq. (4.359) into the last condition in Eq. (4.358) then gives 


^2 -^2 

+ a w sx + qw 31 - 2 o^ci = 0 

3x5 3x1 


(4.360) 


However, since ¥ 31 does not depend on <ji , Eq. (4.360) is satisfied if simul- 
taneously 


Cten = 0 


6 2 ¥ 31 + 3 g ¥ 31 

6x^ 3x1 


+ a¥ 31 = 0 


(4.361) 

(4.362) 
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Since (X is not zero, c x must "be zero. 


Therefore, the final form of the characteristic function is 


W = c 2 p x + c 4 p 2 + c 5 jz( + ¥ 31 (x 1 ,x 2 ) (4.363) 


where ¥ 31 is any function which satisfy Eq. ( 4 . 362 ), 

To solve for the absolute invariants, it is necessary to solve 


or, using Eqs. 


dx-, _ dx 2 


£1 £2 

(4. 323 a, b) and ( 4 . 363 ), 


M 


dx t _ dx 2 _ d6 
c 2 c 4 -(c 5 ^+¥ 31 ) 


(4.364) 


(4.365) 


Eq.. 


Let us next consider the case where W 31 is a function of x x alone, then 
(4*362 ) is reduced to 


dxf 


+ a¥ 31 


0 . 


(4.366) 


The absolute invariants are the two independent solutions of Eq. ( 4 . 365 ) 
which are 


J] = x 2 - ZjL x x 
c 2 


(4.367a) 


f(n) 



¥ 3x e 


c 2 - 1 - 


dxi 


(4.367b) 


Putting these transformations into the Helmholtz equation and making use 
of the condition, Eq. ( 4 . 366 ), we get 



+ 2 f * 
C2 




f 


= 0 


(4.368) 
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which is an ordinary differential equation. 


The above is only a special case of all the possible groups of transforma- 
tion which will reduce the original partial differential equation to an ordi- 
nary differential equation. Other solutions to W and the similarity solution 
for those groups can be obtained in a similar manner. We therefore conclude 
that similarity solutions do exist for Helmholtz equation in rectangular co- 
ordinates, in addition to its separability. 

One final remark concerning the boundary conditions is necessary. When 
it is stated that similarity solutions exist, we mean that not only the partial 
differential equation is reduced to an ordinary differential equation, but 
also that boundary conditions can be transformed satisfactorily. For separa- 
tion of variables in the usual sense (see the brief summary earlier), the 
boundary conditions can always be transformed. This is not true for similarity 
transformation. In the example just treated, Eqs. (4. 567 ) and ( 4 . 368 ), the 
only possible similarity solution is when the boundary conditions are given 
as 


0 

II 

at 

*2 - If x i 

0 

II 

at 



where k x and k 2 are constants. By putting W 31 to zero, Eq, (4.367) gives the 
following boundary conditions for Eq. ( 4 . 368 ), 


f(k x ) = f(ka) = 0 


(4.369) 


Conditions of Similarity for a Given Group of Transformation 

As an illustration, consider the spiral group of transformation where the 
spiral group of transformation where the characteristic function is given hy 


W(xi,x 2 ,jZ$,Pi,p 2 ) 

= c lPl + C 2 X 2 P 2 - C 3 (zS , (4.370) 

The conditions for the existence of similarity solution, Eq. (4.346) then 
become 
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2Cn 


.2 c)hg 

dx x 


hihg 




\ + 202X2 { hlhi Sr 


h| 


a dhj 


a Xs 


+ 231x11202 = 


and 


where 


G 0 + GxPi + G 2 P 2 + G3j^ = 0 


Gi 


G, 

,2 dhn dh; 


, _ f,,Z Bh2 >,3 cfch ; 

r G -Cx [hjh £ ^ - h 2 ^ 


, it it + *-*■ (t) % h ( 5 $ 


+ C2X2 


— J. J 

, 2 dh n Sh2 . , dh.2 dho . , , 2 3 2 b .2 

L h2 a£ SI + 2hlhs hi s£ + hlh2 

- 5hi ^ - hi 5%- • 

OX^ uX^l ox^pox^ 

+ = 3 (h lhS ^ . hi |x) -^(^^-hi^) 

(x) L(3h!h 2 |b. + hi|^') + ' : ^r5 h ^r it1 '? ^ 

y oXj_ \ 0X2 OX' 

- 2 c 2 hfh 2 + c 3 hxhal 


0 

(4.571) 

(4.572) 

(^•373) 


(4.3T 1 *-) 
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G 2 = c 


/ov, in dh T Shi , u2 <dh 2 Shn . in 2 in 6 2 h x 

L 2h x h 2 — L -—L. + h x — -s- — -L + h x h 2 — ■ - j - 

\ OX^ OX2 OXj 0X2 OX^OX2 

_ ZV , 2 ^2 _ h 3 ^hg ^ 

1 Sx x Sx 2 1 Sx 1 Sx 2 y 


_ 2 ,_ S 2 hi 


+ c 2 x £ (2h 1 h 2 + hf ^S. ^JL + hlh 2 s-ax 

\ Sx 2 Sx 2 1 Sx 2 Sxs s> 

Xh 2 ^ h 2 _ v, 3 ^hgA 

3hl SX2 ST 1 537 


+ (c 3 -c 2 )(hfh 2 


Sxs 


hf 5^ 


-i-fhfhB Sil 
h x h 2 \ Sx 2 


(x) 


e/ 5h f h2 + h 


3 Sh2 \ 

1 ^ x l/ 


+ c 2 x 2 ( 3h x h 2 


Shi 


bxg 1 S^s J 


hf 


Sh g \ 


- 2c 2 h.fh 2 + c 3 h?h 2 [ 


(4.375 ) 


G 3 = 


a Cl (3hfhI + 3h?hi |^) + ac 2 x 2 ^fofhf ||i. 


+ 3hfh| + ac 3 h.fh| 

°W 


l^ + h ? Sx x 


2c 2 h x h 2 + c 3 hfh 2 


+ e 2X2 (jhflte + hf 

V 0X2 0X2. 


(4.376) 


It should "be noted that for ¥ in the form given in Eq.. (4.370), the condi- 
tions H x = 0, H 2 = 0 and H 3 = 0 are satisfied identically. From Eq. (4.372) 
we conclude that the G 1 s should all "be zero. Thus, we get 


G 0 = G x = G 2 = G 3 = 0 . 


(4.377) 


Equations (4.371) arid (4.377) are the conditions to be satisfied for the h x 
and h 2 for similarity solutions to exist for the spiral group. As long as 
the functions h x and h 2 for a coordinate system satisfy these conditions, 
then the spiral group of similarity transformation exists for that coordinate 
system. 
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As an application of these conditions, Eqs. (4.371) and (4.377)/ let us 
ask if the spiral group of transformation exists for the rectangular coordi- 
nate frame. For such a coordinate system/ "both h x and h 2 are unity. Equa- 
tion (4.371) becomes 


c 2 = 0 . 


(4.378) 


Equation (4.377), with c 2 = 0, are satisfied identically. 

Thus, we conclude that the spiral group of transformation exists for the 
rectangular coordinate frame, since Its h x and h^ satisfys the conditions 
given in Eqs. (4.371) and (4.377) If c 2 = 0. 

The characteristic function is therefore 


w = c±Px - C 3 ^ 


(4.379) 


and the absolute invariants can be obtained by solving 


dx n 


Cl 



and 


x^ = constant 


(4.380) 


(4.381) 


which gives 


T] = x^ and f (r} ) 



and the Helmholtz equation is reduced to 


f" + (a + )f = 0 . 


By substituting into Eqs. (4.371) and (4.377) "the hi and h 2 functions 
for a cylindrical coordinate (h x =1, h 2 = x x ), it is seen that similarity 
solutions do not exist since they do not satisfy these conditions. Other co- 
ordinate frames can be tested in a similar manner. 
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4.5 CONCLUDING REMARKS 


The method given in this chapter can he summarized as follows: Consider 

a partial differential equation 


where 


E(zj_j ■ • • j Zp ) — 0 


(4.382) 


Zi = Xi 
Z2 = 


5 p-i 


-sk 

A x . 


( dx m _ •[_ ) 

A*y 


(Sx, 


-m' 


This equation is said to he invariant under the infinitesimal contact trans- 
formation 


z_^ — Zi + , i 1, . . ° , p 


if the following condition is satisfied 


UF = ii|^-+ : ...+i ^ = o. (4.383) 

OZj OZp 

Since the functions, {j^, in the transformation are expressed in terms of a 
characteristic function, W, Eq.. ( 4 . 383 ) is used to predict the form of W. 

The invariants can then he obtained hy solving the following system of equa- 
tions : 


..... . (4.384) 
?1 ip 

Finally, using the theorems given in article 4.1.5, the number of variables 
can be reduced by one using the invariants as new dependent and independent 
variables . 
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For simultaneous differential equations , the functions in the infinitesimal 
contact transformation are expressed in terms of characteristic functions-, 
where i = l,...,m and m is the number of dependent variables,, 

The present method is seen to be a systematic way of searching for all 
possible groups of transformation which will reduce the number of variables 
by one. For reducing more variables, the same steps have to be repeated. 
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